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1. Introduction 

The concept of Supersymmetric Quantum Mechanics (SUSY QM) embodies an 
algebraic form of transformations of (complete or partial) spectral equivalence between 
different dynamical systems and in this sense it gives the algebraic tools for spectral 
design of new quantum systems from a given set with controllable energy spectra. At 
present, there is a number of reviews [I] [32] devoted to development and various 
applications, mostly of the linear SUSY QM. The very isospectral transformations 
realizing SUSY represent the Darboux (Darboux-Moutard-Crum) transformations [19]- 
|27j . which are known in the theory of ordinary differential equations for a long time. 

The Supersymmetric Quantum Mechanics in its form was formulated for one- 
dimensional space by E.Witten [28] about thirty years ago(see also [29] - [31]). This 
approach has two predecessors: Factorization Method [32] in Quantum Mechanics and 
Darboux transformation [B5]-[25] for the Sturm-Liouville equation in Mathematics. 

The first one was proposed by E.Schrodinger in 1940-ties |32j (see the excellent 
detailed review [33]), and its basic relation is the representation of Hamiltonian in a 
factorized form, 

h(°\x) = -d 2 + U (0) (x) = q + q~, (1) 
where operators q ± are, 

g ± = T d + dW(x); q+ = (q-y d = ^. (2) 

The factorization (TjQ) is in general possible for an arbitrary potential V(x) (see details 
below). The main idea of Factorization Method is to consider together with the 
Hamiltonian (Op) the partner Hamiltonian, also of Schrodinger form, 

(x) = -d 2 + V {1) (x) = q~q + , (3) 

with the same operators q ± . By construction, these Hamiltonians satisfy the intertwining 
relations, 

h®q + = q+h®) q-h® = /i ( V, (4) 

which play the key role in the method under discussion. Just the intertwining relations 
fl4]) allow to obtain the simple relations between bound state wave functions \l/i ' ) and 
of and correspondingly. Indeed, according to (pEJ) an arbitrary wave function 
gives the wave function of and vice versa, 

g -*(°)oc*«; g+tfWKtfjW (5) 

The energy eigenvalues for ^i ^ and ^i 1 are the same, = E { n \ There are three 
options [M] for the spectra of and : 

- spectra coincide completely if neither q + , no q~ have normalizable zero modes; 

- spectra almost coincide: the spectrum of has one bound state less, if the 
normalizable zero mode of q + exists; 

- spectra almost coincide: the spectrum of has one bound state less, if the 
normalizable zero mode of q~ exists. 
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This is an essence of the Spectral Design method based on the super-Hamiltonian 
H with two non-singular real potentials. The supersymmetry results in isospectral 
relations between the super-Hamiltonian components and whereas two non- 
singular real potentials are different. Each of the operators may have one zero mode 
exp (±W(x)), and their normalizability properties depend on the behavior of function 
W(x): the Hamiltonians and are [31] [M| either isospectral or almost isospectral 
(up to one bound state). 

The connection between a given potential V^(x) and a function W(x) provides 
a basis for the link of Factorization Method with Darboux transformation [19] and 
represents the analytic tool of Spectral Design. The substitution of (j2J) into (JTJ) , (j3J) 
results in: 

V<®(x) = (dW(x)) 2 - d 2 W(x); (6) 
VW(x) = (dW(x)) 2 + d 2 W(x). (7) 

In terms of the theory of ordinary second-order differential Sturm-Liouville equations 
(adapted to the Schrodinger form), 

- *<®"(x) + V^(x)¥°\x) = E^¥°\x), (8) 

the Darboux transformation can be formulated as follows. If one knows its particular 
solution ip(°\x) = exp (— W(x)) for the value E^ = 0, the Eq.(jSJ) is satisfied, and 
an arbitrary solution ty(°\x; E^) of (jSJ) can be transformed (Darboux transformation) 
into a solution of the transformed equation of the same form and with the same value 
of parameter E^°\ 

- ¥ l) "(x) + V^ix^ix) = £ (( % (1) (x). (9) 
The initial and transformed solutions are connected, 

$ (1) (x;£ (0) ) oc (d + (dW(x)))^ \x-E^); (10) 

¥°\x;E®) oc (-d + idWix^W&E^). 

It is clear that one obtains a one-to-one correspondence between solutions of (jSJ) and 
( |T0|) . but up to zero modes of operators (=F<9 + (dW)). 

The nonlinear first-order differential equation (Riccati[l2] equation) can be 
solved analytically for a restricted class of potentials V^°\x) only. The class of 
analytically solvable Riccati equations corresponds to the so called exactly solvable 
models in one-dimensional Quantum Mechanics, and latter ones are related to the so 
called shape-invariant [13] potentials. 

The above construction can be easily presented in a supersymmetric form. Indeed, 
the Hamiltonians h^°\ and intertwining operators q ± can be considered as elements 
of matrix 2x2 operators - super-Hamiltonian and supercharges, 
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which obey the (anti) commutation relations - the simplest realization of supersymmetric 
algebra, 

[tf, Q±] = 0; {Q+, Q-} = H; (Q ± ) 2 = 0. (12) 

The relations (fl2|) of the superalgebra express in the compact form the factorization 
([I]) and the intertwining relations (Hj), while the nilpotency of supercharges follows 
directly from their matrix structure. By analogy with the Supersymmetric Quantum 
Field Theory, the component h^°' of H is called the bosonic component, and h^ 1 ' 
fermionic one. Correspondingly, their eigenstates form bosonic and fermionic sectors, 
although this terminology has no physical sense in the present non-relativistic context. 
Nevertheless, the nilpotency of Q ± can be associated with the fermion property - an 
analog of the Pauli principle. 

Thus in the simplest cases (linear SUSY algebras), intertwining of two differential 
operators (for instance, Hamiltonians of one-dimensional quantum systems) by means of 
Darboux operators of first order in derivatives entails their factorization into differential 
multipliers which are formed by the same Darboux operators. However, in general, 
both interrelation between pairs of dynamical operators ("Hamiltonians") with (almost) 
equivalent spectra and structure of operators which generate the spectral equivalence 
are not that simple [H] - [62]. In particular, the closing of a Nonlinear SUSY algebra 
often produces a hidden symmetry generator which helps to partial solvability of 
spectral problem for Schrodinger Hamiltonians. This interrelation in the one- and two- 
dimensional QM is within the scope of the present review. In particular, we survey well 
motivated solutions to the following problems: 

• what is the correct self-consistent way to generalize the main relations of one- 
dimensional SUSY Quantum Mechanics, including factorization and intertwining 
relations, to the case of multidimensional space? 

• in what cases the higher-order Darboux-Crum transformations can be constructed 
with the help of a sequence of intertwining transformations of lower order SUSY 
which relate a chain of (almost) isospectral intermediate Hamiltonians with real 
nonsingular potentials? 

• what are elementary blocks for a nonsingular factorization of intertwining 
operators? 

• in what way the irreducibility of elementary blocks of isospectral transformations 
reveals itself in the SUSY algebra and in the structure of kernels of those 
transformations? 

• are there more complicated SUSY algebras in one dimensional QM which accept 
hidden symmetry generators ("central charges")? 

• what is the relationship between shape-invariance and spectrum generating algebras 
for SUSY QM systems? 

• how SUSY for non-stationary Schrodinger operators produces hidden symmetries 
and spectrum generating algebras? 
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• is it possible to find pairs of (almost) isospectral scalar two-dimensional 
Hamiltonians intertwined by higher order supercharges, i.e. to avoid the matrix 
components of super-Hamiltonian by means of using Nonlinear SUSY algebra? 

• does the Nonlinear two-dimensional SUSY Quantum Mechanics give the 
constructive way to obtain new exactly solvable and/or quasi-exactly solvable 
models not amenable to the conventional separation of variables? 

The structure of the paper is as follows. After a short reminder of notation and 
basic definitions of SUSY theory of isospectral transformations given above in the 
Introduction, the generalizations of linear SUSY QM for arbitrary space dimensions [34J , 
[64] are considered in Section 2. All of them can be called as direct generalizations since 
they are based on the supercharge operators of first order in derivatives. In particular, 
Subsection 2.1 includes a detailed description of two-dimensional generalization [34] in 
the form similar to the Factorization Method of E.Schrodinger [32], [S3], and Subsection 
2.2 contains two illustrative examples: Pauli equation for fermion motion on the plane 
with external magnetic field [65] - [67] and supersymmetric extension of two-dimensional 
attractive Coulomb potential [34]. In Subsection 2.3, the well developed formalism 
of Supersymmetric Quantum Field Theory is used to present the direct (first-order) 
generalization [64J of Witten's one- dimensional SUSY Quantum Mechanics onto the 
space with arbitrary number d of dimensions. The main element of construction is 
the d— component real scalar superfield depending both on time and pair of Grassmann 
variables. This multidimensional generalization allows also the alternative interpretation 
as a model of d interacting particles on a line [18J. In Subsection 2.4 the structure 
of obtained super-Hamiltonian is analyzed [64] in the Hilbert space which is a direct 
sum of subspaces with fixed fermion numbers. In this basis, the super-Hamiltonian 
has the block-diagonal form with matrix blocks of different matrix dimensionality. 
The supersymmetric intertwining relations between different matrix components of the 
super-Hamiltonian are derived providing their isospectrality. Subsection 2.5 gives a 
brief description of the illustrative example [68] : d = 3 SUSY Quantum Mechanics 
with specific Yukawa-like superpotential provides relations between different sectors of 
nucleon-nucleon and nucleon-antinucleon interactions. 

Section 3 is devoted to building of the Polynomial SUSY algebra for d — 1 with the 
help of dressing chain or ladder construction [44]. It is elucidated how they entail the 
SUSY algebra representation by differential intertwining operators of higher-order in 
derivatives. For algebraic formulation its key ingredient is given by SUSY Hamiltonian 
projections on the formal zero-mode space (the kernel) of supercharge operators which 
are finite-dimensional constant matrices [621 [69] . However the basic elements of spectral 
design in Polynomial SUSY consist not only of intermediate isospectral transformations 
of first order in derivative (of linear SUSY's) but also of the transformations of second 
order in derivative [45j [52} [70] . Three irreducible classes of binary SUSY algebras are 
revealed and it is shown how they serve to complete the tool kit for spectral design. 
In particular, they may help to add or remove one or two adjacent levels for excited 
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states. Thus the classification of irreducible (almost) isospectral transformations and 
related SUSY algebras is outlined (described first in [T4] and systematically in [7T| 172"]). 
In connection to this classification the important problem of a possible redundancy in 
supercharges which can be eliminated without any changes in the super-Hamiltonian is 
analyzed in Subsection 3.4. With the help of a (redundant) multiplication of supercharge 
by a polynomial of super-Hamiltonian it is elucidated that certain irreducible second- 
order supercharges can be embedded into a higher-order SUSY which in turn is 
decomposable in a chain of linear SUSY's. In Section 4 we formulate basic theorems 
which give exhaustive answers on the structure of polynomial one-dimensional SUSY 
algebras and on minimization of this algebra up to its essential part (proofs of these 
theorems can be found in [69 | I7TI 172] ). Apart from the above mentioned trivial ways to 
generate several SUSY algebras for the same super-Hamiltonian there exist possibilities 
to obtain two non-minimizable SUSY algebras with supercharges of even and odd order 
in derivatives (69J [73] . Their combination forms a non-trivial hidden-symmetry operator 
which cannot be expressed as a polynomial of the super-Hamiltonian. Such systems and 
algebras are described in Section 5. They are governed by (reflectionless) potentials of a 
special form arising also as soliton solutions of some non-linear equations [26j E7J . The 
related SUSY algebra can be qualified as Af = 2 superalgebra with central charge made 
of hidden symmetry operators. 

In Sect. 6 the concept of shape invariance is formulated for d — 1 polynomial SUSY 
algebras with the help of intertwining with shift |73j. The arising spectrum generating 
algebras represent intrinsic features of shape invariance and they can be considered as 
deformations of hidden symmetries studied in Sect. 5 . It is elucidated on and exemplified 
for the second- and third-order SUSY shape invariance. The relations with deformed 
Heisenberg algebras [TH [56] are described as well. In the case of third-order SUSY its 
remarkable connection to the Painleve-IV equations [74] is demonstrated (Sect. 6.4) 
and the typical patterns of energy spectra are listed. 

Section 7 is devoted to extension of spectral design methods on nonstationary 
Schrodinger equations [75J. We follow the scheme of SUSY intertwining relations and 
show that addition of time dependence leads to certain restrictions on potentials which 
are assembled in a super-Hamiltonian. In all cases the closure of the SUSY algebra 
reveals a hidden symmetry which becomes a typical rather than exceptional feature 
for non-stationary Hamiltonians. For linear SUSY (Sect. 7.1)this hidden symmetry 
entails simply full separation of variables. For the second-order SUSY (Sect. 7.2) one 
discovers the spectrum generating algebra previously found for the third-order shape 
invariance and this algebra is completely encoded in time-dependent hidden symmetry 
operators. Thus hidden symmetries for nonstationary Schrodinger super-Hamiltonians 
give the universal framework for both stationary hidden symmetries and for spectrum 
generating algebras of shape invariant SUSY systems. 

In Section 8, the idea of one- dimensional Polynomial SUSY QM is generalized 
onto two-dimensional case providing isospectrality of scalar Schrodinger Hamiltonians 
without any intermediate matrix components [76] - [93]. Subsection 8.1 gives the basic 
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formulas [76] - [79J of this approach for arbitrary form of highest (second-order) terms 
in supercharge. The most promising case of Lorentz (hyperbolic) form of second-order 
terms is studied in Subsection 8.2 with special attention to reducible (with twist) and 
irreducible supercharges. The fundamental difference of two-dimensional Polynomial 
SUSY QM is that both SUSY-partner Hamiltonians are completely integrable: they 
commute with symmetry operators which are of fourth order in momenta (Subsection 
8.3). Section 9 is devoted to two new methods [T3], [T6] of solution of two-dimensional 
Schrodinger equations which are not amenable to conventional method of separation 
of variables. Both new methods are based on the isospectrality of scalar Hamiltonians 
established in the previous Section and on the property of shape invariance (Section 6). 
These methods can be considered as supersymmetric separation of variables. The first 
of them (Subsection 9.1) provides the quasi-exact solvability of several two-dimensional 
models, i.e. it allows to find analytically a part of spectrum and corresponding wave 
functions. The second one (Subsection 9.2) gives the exact analytical solution of the 
models but it works only for particular values of parameters of the models. Both 
methods can be applied to several potentials not amenable to standard separation of 
variables, and here they are illustrated by solution of two-dimensional generalization 
of Morse model. In Section 10, we list some related problems of SUSY QM which are 
outside the scope of the present review, and we mention some related open problems 
and perspectives for future study. 

2. SUSY QM in arbitrary space dimension 

2.1. Two-dimensional SUSY QM with first-order supercharges 

The generalization of one-dimensional SUSY QM onto higher space dimensionality 
d > 2 in the Schrodinger equation is of great interest both from the conceptual and 
practical point of view. The forthright attempt to preserve the one-dimensional formulas 
given above but with superpotentials depending on multidimensional coordinate vector 
x = (x%, X2, Xd) leads to potentials which are amenable to standard separation of 
variables (see details in [91]). In what follows such generalizations will not be considered 
since the corresponding Schrodinger equation is reduced to the simpler problems of lower 
dimensionalities. 

A more advanced generalization of one-dimensional SUSY QM onto multidimen- 
sional case was made in the framework of Factorization Method in [31] (see also [95] and 
|96j). The main idea is in involving vector intertwining operators qi(x); I = l,2,...,d, 
which are still differential operators of first order in derivatives, 



In terms of these operators, the initial Hamiltonian h^- ' (x) can be quasifactorized, 



q±(x) = Tdi + (diW(x)); x = (x u x 2 , ...,x d ); d t 



_d_ 

dx{ 



(13) 



(£) = - A (2) + U (0) (x) = g+gf; A (2) = dA 



1, 



(14) 
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where the sum over repeated index I is implied. In the framework of this generalization, 
which can be named as "direct generalization", the Schrodinger operator is related 
to the chain of d other Hamiltonians of Schrodinger form, but (d — 1) of them with 
matrix potentials (see Subsection 2.5 below). 

Now we shall describe this construction [34] in the simplest - two-dimensional d = 2 
- case, where the chain of Hamiltonians includes three species, the original and two 
others, and h^ 2 \ both of the Schodinger form as well. The first of them has 2x2 
matrix potential Vf£\x), while the second one has a scalar potential V^\x). These 
partner Hamiltonians (superpartners) are expressible in a quasifactorized form as well, 

h ( 2 = qiqt + PTPt = -<M 2 + ((diW(x)) 2 - dfW(x)^j + 2d i d k W(x); (15) 

= pfp~ = -df + V^ 2 \x) = -df + (dtWix)) 2 + dfW(x); (16) 

The operators qf were defined in f|T3|) . and the new vector operators pf are also of first 
order in derivatives, 

pf = timl, ( 17 ) 

where is a standard antisymmetric tensor with = +1. 

The chain of two-dimensional Hamiltonians h^°\ h^J , is analogous to the 
one- dimensional partnership h^°\ h^ l \ since its components are also related by the 
intertwining relations with operators q l , pf playing the role of intertwining operators, 

^2W=pr^ (2) ; PthSt = h<?>pt. (19) 

By definition, the intertwining operators q l , pf, and, correspondingly, two terms in the 
definition f lT5|) of are mutually orthogonal, qfpf = 0. These two terms in were 
necessary to obtain the Hamiltonian of the Schrodinger form, and the orthogonality is 
important to provide the intertwining relations ( TT8|) . (170]) . 

The energy spectra of and are, in general, non-overlapping, but the 
intertwining relations f lTHj) . f fT9|) provide the equivalence of energy spectra between 
a pair of two scalar Hamiltonians and 2x2 matrix Hamiltonian hfj^ . The 

"equivalence" means coincidence of their spectra up to zero modes of operators qf,pf. 
Thus, the supersymmetry (supersymmetric transformation) allows to reduce the solution 
of spectral problem for the matrix Hamiltonian h^f to solutions of the couple of scalar 
spectral problems h^°\ realizing the SUSY diagonalization method. Due to the 
intertwining relations, the vector wave functions of matrix Hamiltonian are also 
connected (up to a normalization factor) with the scalar wave functions of scalar 
Hamiltonians h^ \h^ 2 \ 

*f \x; E) = qr¥°\x; E); % = 1, 2; ¥°\x; E) = qf^ ( f\x; E) 
^f\x- E) = pr¥ 2 \x; E); ¥ 2) (x; E) = pf^ 1] (x; E). (20) 

The two-dimensional generalization ( 1T41) . f lT5l) . (1IB]) . ( |T8l) . ( TT9l) of Factorization 
Method can be also formulated compactly [64] in terms of exactly the same superalgebra 
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(|12p . Indeed, three Hamiltonians h^°\ , are combined in the super-Hamiltonian 

in the supercharges Q ± , correspondingly, 

/ o\ 



if, and intertwining operators q l , 



H 



( m 



V 



h 





(i) 
)fc 










\ 



J 



Q + = (Q- 



V 



o 
o 

pf 



o 



pi 








(21) 



-fermionic 



all matrices are 4x4. The components are bosonic, and h\ 

component of the super-Hamiltonian. 

A simple rearrangement of this representation gives an idea of SUSY diagonalization 
of matrix Hamiltonians in terms of the conventional SUSY algebra ( fill . 

U°) \ no) _ ( ^ (0) 

o h,w r o 



H 



h 



Q + = (Q- 




q 



pi pi 



(22) 



2.2. Examples 



The Schrodinger operators with matrix potentials are not something very exotic in 
Quantum Mechanics. In particular, the described two-dimensional generalization of 
SUSY QM was successfully used [65] - [67] to investigate the spectra of stationary Pauli 
operator for fermion in external nonhomogeneous electromagnetic fields. 

1. The Pauli Hamiltonian for nonrelativistic spin 1/2 particle in three-dimensional 
space is written as follows, 



Hp = {id a + eA a f - fia a B a + U, 



(23) 



where a = 1,2,3, e and /x are the electric charge and magnetic momentum of the 
particle, {A a (x)} is the vector part of electromagnetic potential, B{x) = rotA(x) is the 
corresponding magnetic field (in general, nonhomogeneous), U(x) is the scalar potential 
(which may include the electrostatic term eAo(x)), and a a are the standard Pauli 
matrices. Restricting the problem to the case of external potentials which do not depend 
on one of coordinates (on X2 below), the Pauli Hamiltonian ( 1231) . being 2x2 matrix 
operator, can be identified with the matrix component h^; (i,k = 1,2) of the super- 
Hamiltonian in the previous Subsection. Indeed, the fermionic (two-component) wave 
function has to be factorized into the trivial plane wave along X2 with some momentum 
k and nontrivial two-dimensional wave function, 



typ(x) = exp (ikx2)^(x\, x^). 
In the subspace of such functions, Pauli operator can be rewritten as, 



(24) 



H P ( Xl , x 3 ) = -(di+ ieAxf + U + (k + eA 2 f - fia a B a ; 1 = 1,3; A x = A^xx, x 3 ).(25) 
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Comparing fl25|) and f lT5|) . one concludes that their identification (up to replacing 
(xi,x 3 ) 4-> (xi,x 2 )) is possible if external fields in (J2"5]) are expressed in terms of the 
only function W(x\,x 3 ), 

-fiB= (2d 1 d 3 W(x 1 , x 3 ), 0,(dl- dl)W( Xl ,x 3 )). (26) 
Let for simplicity the external sources of magnetic field are absent, 

rotB = divB = (27) 
Then, the function W is a polynomial of fourth order, and magnetic field is of the form, 

— fiB = ((2axiX 3 + cxi + bx 3 + g), 0, (ax\ — ax\ + bx\ — cx 3 + d)), (28) 

with arbitrary constant coefficients a,b, .... The scalar potential U(x\, x 3 ) is a polynomial 
of sixth order in coordinates (see details in [55]) - [57]. Due to results of the previous 
Subsection, the spectral problem for matrix2x2 Pauli Hamiltonian is reduced to a couple 
of spectral problems with scalar Hamiltonians h^°\ h^ 2 \ which is much simpler. This 
example demonstrates the opportunity of dynamical supersymmetric diagonalization of 
a class of matrix problems by means of two-dimensional SUSY Quantum Mechanical 
approach. The identification of Pauli operator Hp with matrix component of super- 
Hamiltonian is possible for much wider class of external fields, if some unitary rotation 
is preliminary done (see details also in [65J - [67J). 

It is instructive to notice that this method works for a wide class of non- 
homogeneous magnetic fields which can be directed at any angle with respect to plane, 
including the field parallel to the plane. The magnetic dipole moment \i of the particle 
is also arbitrary. It is necessary to mention the different supersymmetric approach to 
solution of Pauli equation for the fermion in the plane [97], [HE], [99], [lQQj . which 
works only for gyromagnetic ratio equal 2 and only for (non-homogeneous) magnetic 
fields directed orthogonally to the plane. The latter condition in itself provides the 
diagonalizability of the problem, but SUSY allows to study the symmetry properties 
and the spectrum of Pauli operator. The method of [97], [98], [99] basically has the 
form of one- dimensional SUSY Quantum Mechanics but with matrix 2x2 supercharge 
operators Q ± (see also jlOlj ) . 

2. The second illustrative two-dimensional model [M] corresponds to the initial 
potential of Coulomb form, 

V^\x) = --; r 2 = x\ + xj; a > 0. (29) 

In this case, the superpotential can be taken as W = 2ar, and therefore, the component 
has no bound states at all, 

V^(x) = +-. (30) 
r 

These specific circumstances lead to equivalence of the spectra of matrix Hamiltonian 
operator and scalar component of super-Hamiltonian h^(x). Since the discrete spectrum 
of the latter is well known, the spectrum of the matrix Hamiltonian with potential, 

V^(x) = -^(2^-5 ik ). (31) 
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is known as well. Its wave functions can be obtained from the wave functions of scalar 
Coulomb problem by means of intertwining relations ( 1181) with operators q[ . 

2.3. SUSY QM in arbitrary dimensionality of space as (0 + 1) — dimensional Quantum 
Field Theory 

Both one-dimensional Factorization Method and its two-dimensional generalization were 
shown above to realize the simplest superalgebra ffl2|) of Supersymmetric Quantum 
Mechanics. This observation gives the idea how we could generalize this approach further 
for arbitrary dimensionality d of space: it is necessary to build new realizations of the 
same SUSY algebra (|T2|) . For this purpose it is helpful [64] to exploit some methods 
and experience of well-developed Supersymmetric Quantum Field Theory |1U2] . Indeed, 
new such realizations of the SUSY algebra was obtained [64J by study of some simple 
models of SUSY Quantum Field Theory. The simplest QFT, which can be interpreted 
as nonrelativistic Quantum Mechanics, is the theory in (0 + 1)— dimensional space, 
i.e. fields depend on one variable - the time t. Keeping in mind, that supersymmetry 
must be a property of the model, we consider the superfield, which depends both on one 
bosonic coordinate (time t) and a pair of mutually conjugate Grassman variables 9, 9. 
The dimensionality d of SUSY Quantum Mechanics on target is provided by choosing 
the fields real, but with d— components. 

Thus, the main ingredient of the model is [64J the d— dimensional real superfield, 

W(t;9,6); I = 1,2, ...,d; <pi = fr; (32) 
{9, 9} = {9, 9} = {9, 9} = [t, 9] = [t, 9} = 0. (33) 

As it is typical for the genuine superfield theory, the infinitesimal supertransformation 
in the space {t; 9, 9) is defined as, 

t ->• t' = t - i(9e - e0); 9 -> 9' = 9 + e; 9^9' = 9 + e, (34) 

with a pair of small Grassman mutually conjugate parameters of transformation 
e, e. The infinitesimal transformations ( 1341) can be generated by differential operators 
(supergenerators) in the superspace (t;9,9), 

(35 » 

which produce the finite transformations by, 

U = exp[i(eQ + Qe)}. (36) 
Operators Q, Q satisfy the anticommutation relations, 

{Q,Q} = {Q,Q} = 0; {Q,Q} = i^. (37) 

Due to nilpotency of 9, 9, the bosonic superfield can be expanded as, 

Wit; 9, 9) = Xl (t) + i9^{t) - i^(t)9 + 99D t (t), (38) 
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where xi(t),Di(t) are bosonic real d— component fields, and ipi(t), V>i(i) — mutually 
conjugate d— component fermionic (anticommuting) fields. Restricting ourselves to 
superscalar fields, 

V ' l (t' ] 9',9') = ^(t;9,9), (39) 

one obtains the following rules of transformation, 

5xi(t) = ie^(t) - i^(t)e; 5D t (t) = e0,(t) + ^,(t)e; 

5^ (t) = HD l (t) - m (t) ; 5fa (t) = -ieD t (t) - ex, (t) , (40) 

where dot means derivative over t. 

As usual, for construction of superinvariant action functional we need the operators 
of supercovariant derivatives, 

d-d d d 

D = !»-**• D = m- a W (41) 

which allow to build the superinvariant action, 

S = J dtL{(p) = J dt J d0£(-D<pD<p + W(<p)) (42) 

with arbitrary functional W((p) - superpotential of the model. Expanding it around the 
point xi(t), integrating the action ( )42l) over both Grassmann variables and substituting 
nondynamical variables Di, Di from the corresponding equations of motion, one obtains 
the action as functional of component fields and superpotential, 

S = J dt(x lXl + - 4^ - (dWimdiWix)) - {d l d m W{x))[i, l {t) : Vm(*)]). (43) 
The corresponding classical Hamiltonian is, 

H d = PlPl + {djW){djW) + (didnWftMm], (44) 

where pi is the momentum p\ = 5L/5xi = X/. According to Noether theorem (in the 
supersymmetric context), supercharges in the present field theory are |103j . 

Qti = 5 -§ = (m + (WM; Qa = ^ = Hft + (WM = W + ) f - ( 45 ) 

To present point, the model was formulated on the classical level, but now 
everything is prepared for its quantization (the prescriptions for canonical quantization 
of mixed systems both with bosonic and fermionic degrees of freedom see in |104j ). 
After quantization, coordinates xi, momenta pi and fermionic variables become 
the operators xi,pi,b~f ,bj~ , correspondingly, 

[Pi, x m ] = -iSi m ] {b'l, b m } = S lm ; {bl, 6+ } = {6f, b m } = 0, (46) 

and anticommuting operators bf can be interpreted as creation and annihilation 
operators for the system of d localized scalar particles, but with Fermi statistics. The 
quantized super-Hamiltonian has the form, 

H = Mi + + d^Wb+b-, (47) 
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and quantum supercharges are, 

Q ± ={±ip l + (d l W(x)))bf. (48) 

These quantum operators satisfy all relations of superalgebra (112]) . and they provide 
the required generalization of one- and two-dimensional supersymmetric models of 
Introduction and Subsection 2.1 to the case of arbitrary dimensionality d. 



2.4- The structure of super- Hamiltonian in SUSY QM for arbitrary dimensionality of 
space 



In order to understand the structure of d— dimensional quantum super-Hamitlonian (T4j 
and supercharges (148]) . it is useful [64] to use the fermion number operator, 

N F = b+b^, (49) 

which commutes with H, Q ± as follows, 

[H, N F ] = 0; [N F , Q ± ] = ±Q ± . (50) 

If the basis in the Hilbert space % of states is chosen so that the space is a direct sum of 
subspaces T-L n with fixed number n of fermions (n = 0, 1, 2, d), the conservation law 
()50p means that super-Hamiltonian has the block- diagonal form in the fermionic Fock 
space, 

\ 




H 



( 




















h (d) 



(51) 



/ 



with components acting in the subspace H n . The convenient basis in H n is the basis 
of occupation numbers - C% vectors \ni,n 2 , ■■■,n ( i > with n\ = 0, 1, and Yl l i=i n i = n - I* 1 
turn, the supercharge Q + has the block structure with nonvanishing matrix C^ +1 x C% 



blocks q^n+i below the main diagonal, and Q~ 
the main diagonal, 



matrix C n d x C% +1 



or 
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Qd,d-l 





(52) 



The algebra ( fl2|) and matrix representations ( |5T|) . (|52|) provide the following relations, 

= Qn-l^Qn^n-l + < in+l,n ( ln,n+l = k. + h ^ (53) 

(quasifactorization - compare with (IT5l)). 



qn-l,n-2ln,n-l = ^ hPhP = = 0, 



(54) 
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- itPi = 
a~ h {n+1) 

Hn+\,n'-L 



h {n) q. 



n+l,n' 



in Subsection 2.1), 
= 0; 

(55) 



(orthogonality - compare with the orthogonality pfqj 

h (n+l) + _ + L(n) _ i.(n) - 
LL Hn,n+l — Hn,n+l!± ~ Ik Hn+l,r 

Ik Hn,n+l — Hn,n+lLk i H n +l,n'-k 

(intertwining relations analogous to (jlj)). 

Thus, each matrix component hS n \ (n = 1, ...,d — 1) of the super-Hamiltonian H 
contains two mutually orthogonal terms h^, hr n . These terms are intertwined with 
analogous terms of neighbor matrix components, providing their spectral equivalence 
(again, up to zero modes of supercharges). Due to orthogonality (I54|) . each matrix 
component h^ n \ as a whole, is also intertwined with neig hbor h( n -V and h^ n+1 \ As for 
connection of the wave functions, those of are connected, analogously to (jSJ) and 
( 12 Op with wave functions of neighbor Hamiltonians h^ n+1 \ h^ n ~ l \ This scheme can be 
illustrated by the following diagram, 

[h<®\ <=► \hM + <z» [A® + ^)] <=> ...^-D + ^-i)] <-> h®. (56) 

Thereby, the structure of elements of SUSY QM for arbitrary dimensionality of 
space is known. In order to obtain the explicit expressions for potentials, one has to 
choose any matrix representation for the fermionic creation and annihilation operators 
b^. In particular, choosing 



(57) 



the relations of the Introduction for d = 1 are obtained. For the space dimensionality 
d = 2, the suitable representation for fermionic creation and annihilation operators 
± I — 1, 2 is 4 x 4 matrix, 
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(5* 



leading [64J to the same 4x4 matrix dimensionality of the supercharges and the 
super-Hamiltonian H given by (j2"T|) . 



2.5. The particular case of d = 3 : nucleon-nucleon interaction 

For physically interesting case of d — 3, three creation and annihilation operators 
bf, I = 1, 2, 3, and the supercharges are of matrix dimensionality 8x8. The super- 
Hamiltonian H includes [01] four components: two scalars h^°\ and two 3x3 
matrices hf k , hf^ . According to the general scheme of previous Subsection, there are 
intertwining relations in the pairs ^ h\^; ^ hf^ ^ hP\ providing the 
relations between corresponding scalar and vector (three-component) wave functions. 
In this case, the same commutation and anticommutation relations of the superalgebra 
(OGD are fulfilled. 
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The possibility of using of this 8x8 realization of ffl2|) can be illustrated by 
derivation [68J of simple connections between different sectors of nucleon-nucleon 
and nucleon-antinucleon interactions in the OPEP (one-pion-exchange-potential) 
approximation. The general form of the Hamiltonian of interaction of two non- 
relativistic particles of spin 1/2 with scalar, spin-spin and tensor forces is: 

#(l/2<g>l/2) = -iA( 3 )+V s (r) + i(l+a\a 2 )V CT (^^ (59) 

where <?i, cx 2 are Pauli matrices describing spin operators of two nucleons. Separating 
fl59|) onto singlet spin (S = 0) and triplet spin (5=1) parts, 

V s =° = V s (r) - V a (r); V$ =1 = (V s + V a + 2A^V T )5 lk - dd^Vr, l,k = l, 2, 3, (60) 

and choosing suitably the spherically symmetric superpotential W(r), one may identify 
the first part V s=0 with the scalar potential of the component of the super- 
Hamiltonian: 

V s=0 = -\(dW(r)) 2 - A (3) mr)l = -\(W') 2 - W" + -W'l (61) 

2 2 r 

where prime means derivative over r. Simultaneously, it is possible [68J to identify the 
triplet spin part V^f =1 with the 3x3 component of the super-Hamiltonian, i.e. to 
put: 

V T (r) = V(r); V a (r) = l^W^); V s (r) = \{dW{r)f - \^W(r). (62) 
6 3 2 6 

In the OPEP-approximation, the explicit form of NN interaction potential is: 

yNN = f2{m){did){B2d) (63) 

r 

where fx, f 2 are isospin operators of nucleons, f 2 is a coupling constant, and \i— pion 
mass. The expression (1631) corresponds to the specific forms of tensor, spin-spin and 
scalar terms: 

v nn = 1 /2(^ 2 )£^ ; v nn = 2 A®V» N ; V g NN = -A^V™. (64) 

3 r 

Choosing the superpotential W in the super-Hamiltonian H as 

e fJ.r 

W(r) = 2(f 1 f 2 )f — , (65) 
r 

one can identify [HE] the Hamiltonian H NN with (f) with exponential accuracy 
[W') 2 /W" (see [105]). 

It is important that in the OPEP-approximation, two other components and 
of the super-Hamiltonian H have also simple physical sense (see |68J): they describe 
the triplet spin S = 1 and singlet spin 5 = 0, correspondingly, of nucleon-antinucleon 
interaction H N ^ = hty . Thus, just the parts of H NN and H N ^ with definite total 

spin being the components of the same super-Hamiltonian satisfy the supersymmetric 
intertwining relations of the form f )55p . In case, if there were bound NN states, 
their spectra would (almost) coincide, and we could write the relations between the 
wave functions and ^f^ N , similar to (120]) . In the present case, supersymmetrical 
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intertwining relations allow to connect asymptotic behaviour of the wave functions of 
continuous spectra. The non-stationary and stationary non-relativistic scattering theory 
for a supersymmetrical Hamiltonian was developed in |68j , where the simple connections 
between scattering data - partial amplitudes, phase shifts and mixing parameters - of 
the components of super-Hamiltonian were derived. 

3. Polynomial SUSY in d=l 

3.1. Supercharges of higher order in derivatives from ladder construction 

The supersymmetric methods of the previous section was based on the superalgebra 
(j!2p which includes: supersymmetry, i.e. conservation of the supercharges Q , 
(quasi)factorization of the super-Hamiltonian components, and the nilpotency of 
supercharges which is provided by their matrix form. Just the supersymmetry (i.e. 
the intertwining relations) guarantees the (almost exact) isospectrality between the 
components of super-Hamiltonian with non-singular potentials, providing the most 
interesting practical results of the method. This is a reason to consider generalized 
SUSY with supercharges of higher order in derivatives, and therefore, without quasi- 
factorization of Hamiltonians in terms of components of supercharges. In the context 
of Darboux transformation [19] , such a generalization of the intertwining relations (jlj) 
corresponds to the so called Crum transformation |21j : 

hW<& = dhM, q- N h^ = h^q- N , (66) 

where Y-th order differential operators are: 

N 

q^ = Y. w k(*)d k , w% = (Tl) N (67) 

with differentiate coefficient functions. Herein let us take them real and the pairs of 
intertwining operators to be mutually Hermitian conjugated = (qj^Y- 

In the literature there are several synonyms for the higher-order SUSY algebra: 
originally it was named as a Polynomial (or Higher-derivative) one [SI H5], later the 
title of Y-fold SUSY has been suggested [62] and a more general term of Nonlinear 
SUSY has been used [1061 1107] in a certain relationship to a nonlinear SUSY algebra 
arising in the conformal QM [108J. 

One way to achieve the non-linear SUSY is in combining of several linear M = 1 
SUSY QM systems in the fermion number representation which are implemented by 
nilpotent supercharges of first order in derivatives ([2]) built from a real super-potentials 
W(x), 

qf = r ± = =f<9 + x( x )i x( x )=9W(x) (68) 

In the latter case the SUSY algebra is completed (1T21) by the appropriate factorization 
([I]), © of the super-Hamiltonian [32| [33]. 

Let us proceed by recursion: from a simple Darboux transformation to a ladder [35] 
or a dressing chain [109J made of several simple Darboux steps. To produce the required 
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transformation operators the two different linear SUSY systems may be overlapped. 
Namely, consider two super-Hamiltonians Hi, i — 1,2, Eq. (jllj) . respectively two sets 
of supercharges Qf with coefficient functions Xi — dWi and supercharge components 
rf = =f<9 + Xi- Let us match two elements of super-Hamiltonians with a constant shift, 

h? = hP + X; xl+Xi = xl-X , 2 + ^ (69) 

Evidently, the ground state energies for the above Hamiltonians Ei,E 2 satisfy the 
inequalities E\ > A > — E 2 . The constant shift of the super-Hamiltonian H 2 — > H 2 + A 
does not break or change its supersymmetry but modify the superalgebra. 

After matching the chain of intertwining relations (jl]) can be assembled to produce 
the final Hamiltonian h 2 only from the initial one, by means of a second-order 
Darboux transformation, 

/4 0) r+r+ = r+/if ) r+ = rf(h 2 0) + X)r+ = rfr+(h { 2 ] + \); r 2 r^hf = (h { 2 1] + X)r 2 r^ .(70) 

Following the spectral design one concludes that the two-step ladder SUSY dynamics 
contains the redundant information, namely, about the intermediate Hamiltonians 
h ( l ] = hf + A. Let us remove it and define the two (almost) isospectral components, 

h® = hf ] + Ai = r+rr + Ai; h w = h { 2 1] + A 2 = r 2 r+ + A 2 ; (71) 
r\r+ + Ax = r\r 2 + A 2 ; (72) 

for the super-Hamiltonian (TTTT) for which we have employed an energy reference shifted 
by arbitrary Ai j2 , A = A 2 — Ai . Then the intertwining relations ( 166|) are identical to 
Eq. ( 1701) with q 2 = r\r 2 = (q 2 )^ and the supersymmetry [H,Q 2 ] = [H,Q 2 ] = is 
generated by the conserved nilpotent supercharges, 

<#=f °- o)=^) f , {Qt)2 = m * = °" (73) 

Now in virtue of (ITT]) the algebraic closure of the SUSY algebra is given by, 

{Qt, Qn=( rt4 ^ r - r -° ri+r + )=(H- \i)(H - A 2 ), (74) 

generalizing Eq.( TT2|) . Thus we have obtained the second-order Polynomial SUSY algebra 
[44j-[49j as a concise form of isospectral deformation of a potential system realized 
by a ladder [3TJ |^TJ or a dressing chain |109] - [TT2"] of a couple of one-step Darboux 
transformations or, equivalently, by a second-order Crum-Darboux intertwining [2TJ or 
by a blocking of two linear SUSY with partial overlapping of super-Hamiltonians [H] 
("weak SUSY" |113j ). or by a tower of para-SUSY transformations [114] — [TT8] (see how 
to proceed in |14j). 

The formal zero-modes^l of intertwining operators q 2 form the basis of a two- 
dimensional representation of the super-Hamiltonian, 

h^qt<j>t{x) = = qth (1) <Pt(^ * = 1, 2; 

X We mark them as "formal" because they are not supposed to be normalizable, just being solutions 
of the Schrodinger equation but not necessarily belonging to the set of spectrum eigenfunctions. 
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h^ttix) = W(aO; h^fcix) = Stjtjix), (75) 

3=1 3=1 

due to intertwining relations (|66|) . (1701 . In terms of these Hamiltonian projections - 
constant matrices , the SUSY algebra closure takes the polynomial form [69] (see 
also [62], 

{Qt,Q 2 } = det [M - S+] £=H = det [EI - S~] E=H = V 2 {H). (76) 

Thus both matrices have the same set of eigenvalues which for the ladder construction 
( 1741) consists of Ai, A2 . As the formal zero-mode set is not uniquely derived from ( 1751) 
the matrices are not necessarily diagonal. For instance, the equation r 1 f r^0 + (x) = 
has one formal zero-mode <p2 obeying ^^(x) = and another one obeying 
rf4>i = 0; (f>i = r2<f)f{x) 7^ 0. Evidently the zero-mode solution 4>f(x) is determined 
up to an addition of 4>t- When multiplying these linear equations by one proves 
with the help of Eqs. ([7TT). (172]) that 



- \ 2 )<f>i(x) = 0; (hV-\ 1 )4 > i(x) = C4 > t(xy, S-=i A ' , , (77) 




where C is an arbitrary real constant. If Ai 7^ A2 then by the redefinition 
(Ai — A 2 )0^ = (Ai — \2) ( t ) i + C(p2 {%) one arrives to the diagonal matrix S~. However 
in the confluent case, Ai = A2 = A, C 7^ it is impossible to diagonalize and by a proper 
normalization of the zero-mode (ftf one gets the elementary Jordan cell, 



(hW - A)0+(x) = 0; (hU-X)<t>t(x) = <f>i(x); 8~ = * \ . (7 




Therefore in the confluent case the zero-mode <pf is not anymore a solution of 
the Schrodinger equation but it is an adjoint solution which can be obtained by 
differentiation, 0^ = d$\jd\ . Still the intermediate Hamiltonian h = r^rf + A = 
r 2 r 2 + A is well defined and therefore the intermediate isospectral partner <pf(x) of the 
zero-mode 4>\{x) is a solution of Schrodinger equation with the above Hamiltonian. 
The analysis of the matrix S + is similar. Thus we convince ourselves that in 
general the Hamiltonian projection onto the subspace of supercharge zero-modes is not 
diagonalizable but can be always transformed into a canonical Jordan form. 

To complete the description of Polynomial SUSY algebras generated by a second- 
order ladder one should take into consideration also the degenerate case when 
Ai = A2 = A, C = 0. For this choice the matrix S~ is automatically diagonal and both 
zero modes (j)f 2 ( x ) are (independent) solutions of the Schrodinger equation with the 
Hamiltonian h~ . Then it can be proved [69J that the intertwining operator q 2 is just 
a linear function of this Hamiltonian, q% = A — Hence the intertwining is trivial 
fiW — ^(0) anc i sucn supercharges must be ruled out due to their triviality. However for 
SUSY of higher order N > 3 the removal of such blocks from supercharges may lead to 
less-dimensional, ladder-irreducible SUSY algebras (see Subsection 3.3). 
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3.2. Polynomial SUSY QM from a ladder of linear SUSY systems 

Let us give the general description of the iV-step ladder of linear SUSY algebras which 
entails the Polynomial Superalgebra of iVth-order |44j . [4"5] . We introduce a set of 
first-order differential operators for intermediate intertwinings , 

rf = ^d + X i(x), l = l,...,N, (79) 

and the relevant number of intermediate super potentials Xi{ x )- The set of the initial, 
= ho, the final, = hjq and intermediate Hamiltonians, hi = —d 2 + vi(x) consists 
of Schrodinger operators, nonsingular and real ones. They obey the matching (ladder) 
relations, 

hi = r z ~ ■ rf + Xi = rf +l ■ rf +1 + A/+i, I = 1, . . . , N - 1, 

ho = hW = r +. r - + \ i; h N = hW=r„-r% + \ N . (80) 

These ladder relations correspond to the (dressing) chain equations on super-potentials, 

Wi{x), X i(x)=dW l (x), 

Vl (x) = { X i{x)f + ( X i(x))' + \t = {Xl+i(x)) 2 - (xi+i(x))' + A,+i (81) 

The corresponding intertwining (Darboux) transformations hold in each adjacent pair 
of Hamiltonians, 

hi-i • r+ = r+ • hi, rf ■ h t -i = h ■ rf, (82) 
and therefore the chain of iV overlapping SUSY systems is properly built, 

[H h R+] = [H h Rf\ = 0, Hi - Xi = {R+, Rf}, (83) 

Now let us omit a chain of intermediate Hamiltonians between and and 
produce the Higher-derivative ~ Polynomial ~ Nonlinear SUSY algebra for the super- 
Hamiltonian H given in Eq.flTTJ). The intertwining between and W- 1 ' is performed 
by the Crum-Darboux operators, 

q + = r + . . . r+, q~ = r„ . . . rf . (84) 

The SUSY symmetry [H, Q^] = [H, Qjj] = 0, is accomplished by the supercharges 
of the same matrix structure (1731) and the super-Hamiltonian is represented by finite- 
dimensional matrices on the subspaces of supercharge zero-modes, 

q ±<Pt(x) = 0; z = l,2,...,iV; 

N N 

h {0 Ui(x) = £ S^j(x); h^Ufix) = S ^t(x) (85) 

0=1 3=1 

due to intertwining relations (|66|) . In terms of the constant matrices , the algebraic 
closure is given by a non-linear SUSY relation [691 E2] , 

N 

{Qti,Q N } = det [EI - S+] E=H = det [El - S~] E=H = V N (H) = J] (H - X t ) . (86) 

i=i 
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Again both matrices have the same set of eigenvalues which for the ladder construction 
(J8"0j) consists of Ai,...Ajv • If the degenerate roots appear the matrices 5 ,± are not 
diagonal and consist of Jordan cells. If all intermediate hi are hermitian, nonsingular, 
then A; are real and each ladder step is well defined. 

However it turns out that not all higher-order intertwining operators can be 
built by a ladder algorithm based on linear SUSY elements with real and nonsingular 
intermediate Hamiltonians. 

3.3. Generation of Polynomial SUSY QM by ladder: irreducible SUSY blocks of type I, 
II, and III 

Let us elucidate the circumstances which may obstruct the SUSY ladder 
decomposition of a Polynomial SUSY algebra. The class of Polynomial SUSY algebras 
can be extended admitting complex A; and singular hi. In fact, the full variety of 
building elements for non-linear SUSY can be established within the class of intertwining 
(or transformation) operators qf of second-order in derivatives. Eventually one has to 
find the basic set of Crum-Darboux operators (1671) with nonsingular coefficient functions 
which produce a nonsingular potential after intertwining fl66|) with the Hamiltonian 
possessing a smooth initial potential 

When taking a supercharge of second order in derivative with real coefficient 
functions one always can find formal (not necessarily normalizable) zero-modes of the 
intertwining operators and, further on, the 2x2 matrix representation (1751) for 
the super-Hamiltonian components hS Q \hS^ by matrices S ± . The latter matrices are 
real but, in general, not symmetric. Therefore the first obstruction for the ladder 
decomposition may arise when the reduction to a Jordan form does not yield real 
eigenvalues. For instance, if h^(j)~(x) = Aejfc0~(x) then the eigenvalues of S + = Ae 
are imaginary, mutually conjugated ±i\. The possibility of complex pairs of mutually 
conjugated roots in a Polynomial SUSY algebra can be easily read off from its closure 
( 1751) as for supercharges with real coefficients polynomials V2{H) possess real coefficients. 
We call this kind of irreducibility to be of type I. Its elementary block corresponds to 
the polynomial V^iH) = (H + a) 2 +d, d > and its analytical properties have been 
investigated in [45j, [119j . Some examples of related isospectral potentials are described 
in [T20] . 

Next, one has to ensure the positivity of the SUSY algebra relation ( 1761) in 
a particular differential realization of a super-Hamiltonian H with real non-singular 
potentials and the supercharges Qn,Qn (with N = 2 in our case) made of differential 
operators with real coefficients. Let the energy spectrum Ey, j — 0, 1, . . . ; Ej < Ej + i of 
H be discrete, for simplicity. Then, 

V N {E 3 ) = (Q+v^lQ+v^) + (Q-^IQ-^.) > o, (87) 

if the action of supercharges is well defined in the Hilbert space spanned on 
eigenfunctions of a super-Hamiltonian. It can be extended on a continuous energy 
spectrum with usage of wave packets. 
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Thus for non-singular potentials the acceptable disposition of polynomial roots 
related to zero-modes of a supercharge ensures non-negative values of V(E) for each 
energy level of a Hamiltonian. Accordingly , the following positions for polynomial real 
roots are acceptable(for a pair of complex, mutually conjugated roots the positivity is 
obvious) . 

A. Ai < A 2 < E Q or Ai = E ; A 2 = E x . 

The related SUSY algebra is well embedded into a chain or ladder realization. It is 
reducible because one can gradually add/remove Ai and then A 2 without breaking 
the positivity of intermediate SUSY algebra. The coincidence of roots and energies 
correspond to the isospectral transformation with deleting/ inserting energy levels. For 
instance, if Ai = Eq\ A 2 = E\ then two pairs of zero- modes of q 2 can be chosen as 
solutions of Schrodinger equations with Hamiltonians h' ' , . One can implement [35] 
the energy levels E ,E 1 to appear in any of the Hamiltonians h^\h^> but each level 
only once, either in or in 

B. E < Ai < A 2 < Ei or E = X x , A 2 < E x or Ej < Ai < A 2 < E j+1 , 1 < j. 

A pair of real roots is inserted between adjacent energy levels. Such an algebra cannot be 
decomposed into a chain of two linear SUSY with non-singular intermediate potentials 
as the removal of any of roots Ai i2 immediately breaks the positivity in flHTj) . Then the 
intermediate Hamiltonian acquires inevitably a singular potential which usage would 
lead to the loss of isospectrality. The related Darboux transformations had been known 
in 50ties [23] . We call this irreducibility to be of type II. The examples and certain 
theorems are given in [52], [120j . 

C. E < Ai = A 2 < E x or Ej < \% = A 2 < E j+U 1 < j. 

This is a confluent case which seems to be obtained as a limit of the previous one. 
However, one-dimensional QM does not allow degenerate levels and, besides, the matrix 
projection for the corresponding super-Hamiltonian contains a non-trivial Jordan cell. 
Therefore we specify this separate one. This kind of irreducibility is named to 

be of type III. One may find more information on the analytical properties of related 
potentials in [121j . 

One can use these second-order blocks to build a iVth-order Polynomial SUSY 
system. Their general form is again given by Eq. fl86l) if to accept the presence of 
complex conjugated roots A;. 

Meanwhile it seems that a pair of supercharge zero-modes or even a pair of new 
excited energy levels of the super-Hamiltonian can be always inserted by successive 
application of an appropriate ladder construction described in the previous Section using 
first-order intertwining transformations between regular Hamiltonians. But the order of 
a relevant ladder of first-order transformations and respectively of the final Polynomial 
SUSY will be certainly higher than two. We come to the problem of possible relationship 
between first-order reducible and irreducible SUSY algebras having the same super- 
Hamiltonian. 

The related question concerns the degenerate roots. These roots are distributed 
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between different Jordan cells in the matrices S^. The problem is on how many Jordan 
cells may coexist for the same eigenvalue and what is their role in the supercharge 
structure. 



3.4- Minimization of SUSY algebra for a given Hamiltonian and emergence of 
irreducible blocks of type II and III 

Let us reveal a possible redundancy in supercharges which can be eliminated without 
any changes in the super-Hamiltonian. There exists a possibility when the intertwining 
operators and p N for N > N\ are related by a polynomial factor F(x) depending 
on the Hamiltonian, 

g+ = F(hW)p + Nl = p + N F(h^); q~ N = F(h^)p Ni = p N F(h^). (88) 

Obviously in this case the reduction to the second supercharge does not result in any 
modifications of potentials. 

Thus the problem consists in how to factorize a minimal essential part of the 
supercharge and avoid numerous SUSY algebras generated by means of "dressing" flHHj) . 

Let us illustrate how it could work in the following example: 
the matrix S~ for the minimizable intertwining operator k% with Jordan cells of different 
size having the same eigenvalue. It is generated by the operators, 

p ± = T d + x, h {l) = p'p + + A, k^ = -p + p~p + =p + (X - h {1) ). (89) 

Respectively, the basis of formal zero-modes (eigen- and associated functions) for the 
intertwining operator k^(f)j = generates the non-diagonal matrix S~ , 

p-p+(j)+ = 0+ — > = A0+ + 0+; / A 1 \ 

= — ► hP><f>t = A0+; S- = A . (90) 
^ 0, p-p+rt = — ► = A0+; y A y 

Thus the algebraic redundancy in the operator k£ finds its unambiguous track in 
presence of two Jordan cells in the characteristic matrix S~ with the same eigenvalues. 

The supercharge components cannot be factorized in the form (|88l) if the polynomial 
Vn(x) in the SUSY algebra closure (1761) does not reveal the degenerate zeroes. Indeed 
the SUSY algebra closure contains the square of polynomial F(x), for instance, 

k~ N k + N = Fih^p^p+FihU) = F*(hU)V Nl (hV), (91) 

where Vn 1 {x) is a polynomial of lower order, Ni < N — 2. Therefore each zero of the 
polynomial F(x) will produce a double zero in the SUSY algebra provided by (19TI) . 

Thus the absence of double zeroes is sufficient to deal with the SUSY charges 
non-factorizable in the sense of Eq. (1881) . However it is not necessary because the 
degenerate zeroes may well appear in the ladder construction giving new pairs of 
isospectral potentials (see, for instance, [15] for the Polynomial SUSY of second order). 

Now we unravel the origin of irreducible, type-II and -III transformations based 
on the strip-off factorization. Let us consider an example of irreducible SUSY of type 
II with supercharges of second order in derivatives (see previous Subsection). Suppose 
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that it realizes insertion of two new energy levels between the ground and first excited 
states. Then three lowest energy levels E < E\ < E 2 are of importance to study the 
relevant SUSY systems: the ground state level is degenerate between SUSY partners 



E, 



(0) 



Eq^ whereas the two excited levels are present only in the 



and h?>, i.e. 
spectrum of h^ 1 '. 

Certainly one can use the ladder construction fl80l) - fl84|) to prepare the same spectral 
pattern. For this purpose, the intertwining operators (|84p of forth order in derivatives 
must be employed. Indeed, one can prescribe the ladder steps for qf as follows: start 
from a pair of isospectral Hamiltonians with ground state energies £"3; generate the level 
E in the Hamiltonian using the intertwining operators rf, then sequentially create 
in the spectrum of the state with energy E2 < E3 by means of rf, next the energy 
level Ei < E2 using rf and finally the ground state with energy Eq < E\ exploiting 
rf. These elementary steps are reflected in zero-modes of the intertwining operators 
= rfrfrfrf and q~ = r^r^r 2 r^ . Namely the ground state of is a zero-mode of 
[i.e. of q~) whereas the states of h^ 1 ' for E , Ei, E 2 are annihilated by the product 
(i.e. by q + ) according to Eq. ( |85l) . As the ground state energies coincide for 
h,(°\h^ the Hamiltonian projections on the qf zero-mode space - the matrices 
are, in general, not diagonalizable but have one rank-two Jordan cell each. Thus, for 
instance, 

C \ 




r 2 r t r 4 



s~ 



I E 








Ei 








E 1 




s- 



/ E 








E ) 

where a non-zero constant C can be normalized to C 
S~ in (1921) is achieved by means of re-factorization q + ■ 
the generation of ground state for is associated now with 
Polynomial SUSY algebra shows up one degenerate root, 

V 4 (H) = (H- E ) 2 (H - E X )(H - E 2 ). 



C 

Eq 









E 2 








E 1 J 



(92) 



1. The canonical Jordan form 



r{ ri ri ri 



r + r + r + r + go faofc 

Respectively, the 



(93) 



Evidently, this fourth-order algebra cannot be optimized to a lower-order one because 
there is no replication of roots in different Jordan cells of matrices S 1 * 1 (see the "Strip- 
off" theorem on minimization in the next Section). However one may perform fine- 
tuning of Darboux transformation parameters to provide the constant C = in f )92|) . 
This peculiar choice provides two rank-one cells in (1921) with the same eigenvalue Eq. 
The SUSY algebra is still given by Eq. ( 193]) but the intertwining operators reveal a 
redundancy, 



qf = (h® - E ) q+. 



(94) 



By construction the left-hand side of this relation is fully factorizable in elementary 
binomials r^~ with hermitian nonsingular intermediate Hamiltonians. But in the right- 



hand side the operator q 2 



r 3 r 4 



does not admit a further factorization with a 



nonsingular intermediate Hamiltonian because after removal of the redundant factor 
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(h^ — Eq) such a factorization is forbidden by the positivity of the SUSY algebra, 
Eq. (JET}. 

One can easily extrapolate the previous argumentation onto the case of additional 
degeneracy of excited levels E\ = E 2 to analyze the irreducible SUSY of type III. Thus 
we come to the conjectures that: 

a) the factorization fl84"j) of intertwining operators is not unique and there exist 
options to have more reducible ladders and less reducible ones with a larger number of 
singular intermediate Hamiltonians; 

b) irreducible algebras of type II and III can be identified with special cases of fully 
reducible higher-order algebras of ladder type when the Hamiltonian projections 
have an appropriate number of pairs of Jordan cells with coinciding eigenvalues. 

We notice that one may save time in calculations when instead of a more lengthy 
binomial ladder construction the isospectral irreducible transformations of type II or 
III are exploited to embed pairs of energy levels between two excited ones. In the 
next section the results of a more rigorous investigation of the relationship between the 
reducible and irreducible intertwinings are presented. 

4. Rigorous and exhaustive results on structure and (ir)reducibility of 
Polynomial SUSY 

Let us summarize the organization of Polynomial SUSY QM in a more systematic way. 

4-1. Basic theorems on the structure of QM with a nonlinear SUSY 

First the following theorem regulates the very structure of Polynomial SUSY in terms 
of the Hamiltonian projection of the kernels of intertwining operators. 

Theorem 1 . 

Let (jy^{x), n = 1, . . . , N be a basis in kerq^: 

qUl = 0, q~ N = (flft)*. (95) 

Then: 

1) the action of the Hamiltonians h^°\h^ on the functions (fin( x ) ^ s described by 
constant N x N matrices, 

N N 

hPH- = Y, S nm<t>7n\ = £ S~ m 0+; n = l,...,N; (96) 

m=l m=l 

2) the closure of the supersymmetry algebra takes a polynomial form, 

{Q + , Q-} = det[M - S + } E=H = det[M - S~} E=H = V N (H), (97) 

where I is an identity matrix and S ± is the matrix with entries S^ m . 
Corollary 1. The spectra of the matrices S + and S~ are equal. 
Now we describe the minimization ("Strip-off") theorem. 
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A basis in the kernel of the intertwining operator in which the matrix S of this 
operator has a Jordan form is called canonical] elements of a canonical basis are called 
transformation functions. 

The potentials V\{x) and Vzix) of the Hamiltonians and h^ 1 ' are interrelated 
by the equation 

V^(x) = V {l \x) - 2[\nW(x)}" 1 (98) 

where W(x) is the Wronskian of elements of an arbitrary (a canonical as well) basis in 
kerq N . 

The intertwining operators are called minimizable if these operators can be 
presented in the form 

q + N = r(h^)p + M = P + M V(hW), q N = V(hW)p M = p M V{h^\ (99) 

where p^ are operators of order M which intertwine the same Hamiltonians as and 
V(x) is a polynomial of degree (N — M)/2 > 0. Otherwise the intertwining operators 
q^ are named as non-minimizable. 

The following theorem contains necessary and sufficient conditions under which an 
intertwining operator q^ = (q^y is minimizable or not (a proof can be found in [69J). 

("Strip-off") Theorem 2 on minimization of an intertwining operator q~^ 
(and equivalently q^ ) 

An intertwining operator can be presented in the form 

m 

ti=p + MU^- h(i) y ki > ( io °) 
i=i 

where p^ is a non-minimizable operator intertwining the same Hamiltonians as q~^ ( so 
that p\/[h^ = h^p^), if and only if a Jordan form of the matrix S + of the operator q^ 
has m pairs ( and no more ) of Jordan cells with equal eigenvalues A; such that, for the 
l-th pair, Ski is an order of the smallest cell and k\ + Ski is an order of the largest cell . 
In this case, M = N — 2 Y2h=i = Ya=i ki> where the ki, m + 1 ^ I ^ n are orders of 
the remaining unpaired Jordan cells. 

Remark 1 . A Jordan form of the matrix S + of the intertwining operator q~^ cannot 
have more than two cells with the same eigenvalue A; otherwise ker(A — h^) includes 
more than two linearly independent elements. 

Corollary 2. Jordan forms of the matrices S of the operators q~^ and q~^ coincide 
up to permutation of Jordan cells. 

If a Jordan form of the matrix S of an intertwining operator has cells of order higher 
than one, then the corresponding canonical bases contains not only formal solutions of 
the Schrodinger equation but also formal associated functions, which are defined as 
follows [T22] . 

A function ip n ^{x) is called a formal associated function of i-th order of the 
Hamiltonian h for a spectral value A n if 

(h - \ n ) i+1 4> n ,i = 0, and (h - A n ) $ 0. (101) 
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The term "formal" emphasizes that this function is not necessarily normalizable (not 
necessarily belongs to L 2 (M)) and therefore it does not belong to the energy spectrum 
and, accordingly, is not included into the spectral resolution of identity. In particular, 
an associated function tp n ^ of zero order is a formal eigenfunction of h (not necessarily 
a normalizable solution of the homogeneous Schrodinger equation). They appear in the 
subspace of zero-mode solutions for irreducible intertwining operators and play role in 
building S-matrix of the Hamiltonian projection. 

4-2. Classification of really (ir)reducible SUSY transformations 

Assume that the intertwining operators are represented as products of the 
intertwining operators k^_ M and < M < N so that 

Qn = Pm^n-mi In = ^n~mPm'i pXi^m — h°^P~M> Pm^ 1 ^ = ^mPm'i 

k^_ M h { - i) = h M k+_ M , k N _ M h M = h ( - l) k~_ M , h M = -d 2 + v M (x), (102) 

where the coefficients kj^_ M and pj^ as well as the potential vm(x) may be complex 
and/or singular. The Hamiltonian Hm is called intermediate with respect to and 
hP~'. In this case, by Theorem 1, the spectrum of the matrix S of the operator is a 
union of the spectra of the matrices S for the operators k N _ M and p^ 4 . 

The intertwining operators are called reducible if these operators can be 
presented as products of two nonsingular intertwining operators (with real coefficients) 
k%_ M and p%, < M < N so that Eqs. ffT02l are valid and the intermediate 
Hamiltonian Km has a real nonsingular potential. Otherwise q^ are called irreducible. 

Irreducible, non-minimizable, intertwining operators of second order with real 
coefficients are divided into three types: [Tj] . 

An irreducible intertwining operator of I type is a differential intertwining operator 
with real coefficients for which eigenvalues of the matrix 5* have nontrivial imaginary 
parts and are mutually complex conjugate. 

An irreducible intertwining operator of the II type is a differential intertwining 
operator q^ (or equivalently q^ ) of second order with real coefficients such that: 

(1) eigenvalues of the matrix S + of the operator q^ are real and different; 

(2) both elements <pt(x) and (p£( x ) °f a canonical basis of ker q% have zeroes. 

An irreducible intertwining operator of the III type is a differential intertwining 
operator q% (or equivalently q^ ) of second order with real coefficients such that: 

(1) the eigenvalues Ai j2 of the matrix S + of the operator q% are equal, Ai = A 2 ; 

(2) a canonical basis in kerg^ consists of formal eigenf unctions, <^ (x), and 
associated functions, cpf^x), of the Hamiltonian h^ 1 ' which assemble into a Jordan cell, 

/>( Vo = Wo, - AiV£ = tf ; 

(3) (pf Q (x) has at least one root. 

Other types of irreducible non-minimizable intertwining operators of second order 
do not exist. 
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Now let us state (in a concise form) two Theorems which characterize reducibility 
of intertwining operators of any order (proven in (TTj [72]) in an exhaustive way: 

Theorem 3. (on reducibility of "dressed" non-minimizable intertwining 
operators) 

It asserts that for any non-minimizable intertwining operator with real spectrum of the 
matrix S, one can "dress" it, i.e. find and multiply by an appropriate polynomial of the 
Hamiltonian, so that the resulting (certainly minimizable) intertwining operator can be 
factorized into a product of intertwining operators solely of first order with intermediate 
real non-singular Hamiltonians. Thus any Polynomial SUSY algebra with real roots can 
be embedded into a higher order algebra corresponding to a ladder of linear SUSY. 

Theorem 4. (on complete reducibility of non-minimizable intertwining 
operators) 

It asserts that any non-minimizable intertwining operator with arbitrary spectrum of the 
matrix S can be factorized into a product of intertwining operators of first order and 
irreducible second-order intertwining operators of the I, II and III type with intermediate 
real non-singular Hamiltonians. Still this factorization is not necessarily unique and, in 
particular, it may contain more or less blocks of first order. 

These theorems have been formulated at full length and rigorously proven in 
[69] [7TI [72] for the class K of potentials V(x) such that: 

1) V(x) is a real- valued function from C|°; 

2) there exist numbers Rq > and e > (Rq and e depend on V(x)) such that the 
inequality V(x) ^ e takes place for any |x| ^ R ; 

3) the functions 



are bounded for x ^ R and x ^ — Rq, respectively. 

5. Hidden symmetry in One-dimensional SUSY QM 

5.1. Simplest SUSY with Hidden symmetry 

The emergence of hidden symmetries happens to be an intrinsic feature of SUSY 
algebra being related to so called central charges. Specifically in SUSY QM there are 
classes of potentials in super-Hamiltonians which structure is dictated by a hidden 
symmetry. In two (and more) dimensions the variety of such isospectral QM systems is 
rather typical and their spectra and eigenfunctions reveal partial or complete solvability 
accounted for by SUSY separation of variables (see Sections 8 and 9). 

However hidden symmetries exist also in one-dimensional SUSY QM if the 
same super-Hamiltonian conserves two non-trivial (non-minimizable) supercharges [73J. 
Namely the two pairs of conserved supercharges K + ,K~ and P + ,P~ form two SUSY 
algebras for a hermitian Super-Hamiltonian H. 



X 




(103) 



±Ro 
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Let us examine the algebraic structure of the simplest non-linear SUSY with two 
non-minimizable supercharges, 

fc ± = d 2 t 2f(x)d + b{x) =f f'(x); 

p ± = T d + X (x), (104) 

induced by a complex supercharge of second order in derivatives (see the preprint version 
of [73] and [69] ). The supersymmetries (j!20p generated by K + , K~ and P + , P~ prescribe 
that 

d 



y(o),(i) 



6 =r --T f + {tf) + Tf-' (105) 

V (0M " = X 2 =F V, (106) 

where x, / are re& l functions and a, d are real constants. The related superalgebra 
closure for K + , K~ and P + , P~ takes the form, 

{K + ,K-} = (H + a) 2 + d, {P + ,P-} = H. (107) 

The compatibility of two supersymmetries is achieved by solutions of the following 
equations 

X = V\ / 2 + £-^V-t4-« = X 2 = 4/ 2 . (108) 



+ - — 

V 








d 







w-- a ' 



2/ \2fJ if 

Eq. (1108p represents a nonlinear second-order differential equation which solutions are 
parameterized by two integration constants. Therefore the existence of two SUSY 
reduces substantially the class of potentials for which they may appear. Evidently 
Eq. (11081) can be integrated into the first-order one, 

(fT = 4/ 4 + 4a/ 2 + 4G / -d = $ 4 (/), (109) 

where G is a real constant. 

The solutions of this equation are elliptic functions which can be found in the 
implicit form, 

/•/(*) Af 

\ , =±(x-x ), (110) 

J to VMT) 

where / and x are real constants. 

They are nonsingular if: 
a) $4(/) has three different real roots and the double root /3/2 is either the maximal 
one or a minimal one, 

$4(/)=4(/-^) 2 ((/ + ^) 2 -(/3 2 -e)), 0<e</3 2 . (Ill) 

Then there exists a relation between constants a, d, Gq in terms of coefficients 0, e, 

a = e-^<0, G = /3(/3 2 -e), d = /3 2 - e) . (112) 
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The constant fo is taken between the double root and a nearest simple root, 
b) < &4(/) has two different real double roots which corresponds in fillip . (I112p to 
Go = 0, j3 2 = e > 0, a = — e/2, d = — e 2 /4. The constant fo is taken between 
the roots. 

The corresponding potentials U*- -*, V^ 1 ' are well known [1231 1124] to be reflect ionless, 
with one bound state at the energy (/3 2 — e) and with the continuum spectrum starting 
from (3 2 . Respectively the scattering wave function is proportional to exp(ikx) with 
k = - p 2 . 

In the case a) the potentials coincide in their form and differ only by shift in the 
coordinate ( "Darboux Displacement" [125J), 

V^'M = p2 2e (0))(1) = 1 f3-^/e 

ch 2 (^(a;-4 0) ' (1) )) 4 v^ /? + ' 

and in the case b) one of the potentials can be taken constant, 

2 

ch 2 (/3(~-x )), 

For these potentials one can elaborate the relations of extended SUSY algebra. 

The initial algebra is given by the relations (JTUTJ). It must be completed by the 
mixed anti-commutators containing symmetry operators, 

{K + , P-} = {K-, P+}t = B(H) + S(H), (115) 

where the first term is Hermitian and in general (see the next Subsection) a polynomial 
of the super-Hamiltonian and the second one is anti-Hermitian and not a polynomial 
of H. In the example under consideration the first symmetry operator turns out to be 
constant, B(H) = Go after taking into account (11041) and (1109j) . Meanwhile the second 
operator reads, 

3__, A„ 3, 



y(0) = y(l) = 01 J- 1 _ _ j ^ (n4) 



£{H) 



Id 3 - (aI + -V(x))a--V'l 



x 



(116) 



in the notations H = —d 2 I + V(x). By construction the operator S(H) realizes a new 
symmetry for the super-Hamiltonian. Directly from Eq. (11161) one derives that, 

- 8 2 {H) = H[(H + a) 2 + d]-G 2 = (H- E b ) 2 (H - f3 2 ), (117) 

where £& = (3 2 — e is the energy of a bound state. Thus (some of) the zero modes of £(H) 
characterize either bound states or zero-energy states in the continuum. However there 
exist also the non-normalizable, unphysical zero-modes corresponding to E = Eb, j3 2 . We 
remark that in the case b) only the Hamiltonian has a bound state. Accordingly 
the physical zero modes of £(H) may be either degenerate [case a), broken SUSY] or 
non-degenerate [case b), unbroken SUSY]. 

The square root in ( j!17jl can be established unambiguously from the analysis of 
scattering (transmission) coefficients, 

8(H) = t(E b - H)^H-f3 2 . (118) 
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We notice that the symmetry operator (11161) . (11181) is irreducible, i.e. the binomial 
(Ei, — H) cannot be removed. Indeed the elimination of this binomial would lead to an 
essentially nonlocal operator. 

When taking Eq. (II 181) into account one finds the mixed anti-commutators of the 
extended SUSY algebra ( 11151) in a non-polynomial form, 

{K + , P~} = {K~, P + y = G - i(H - E b )y^H~^. (119) 

Thus the "central charge" of the extended SUSY is built of the elements (I107p and 
(11191) and cannot be diagonalized by a unitary rotation with elements polynomial in H. 
Therefore the SUSY algebra is extended in the class of differential operators of finite 
order. 

5.2. General case: one-dim SUSY algebra with Hidden symmetry 

Let us examine the very possibility to have several supercharges for the same super- 
Hamiltonian. We remind that a number of supercharges can be produced trivially 
with the help of multiplication on a polynomial of the Hamiltonian (see Section 4). 
Certainly such supersymmetries are absolutely equivalent for the purposes of spectral 
design and one must get rid of them. It was proved in [69J that the infinite set of possible 
supercharges can be always optimized so that no more than two nontrivial supercharges 
remain. They can be used to generate all other ones by "dressing" with polynomials 
of the Hamiltonians. Thus in one-dimensional (scalar) QM one may have the Af = 1,2 
SUSY only. 

Correspondingly suppose that the super-Hamiltonian H admits two supersymme- 
tries with supercharges K and P^ made of differential intertwining operators of order 
N and Ni respectively, 

[H,K ± ] = [H,P ± ] = 0. (120) 

The second pair of supercharges P^ 1 is assumed to be made of differential operators of 
lower order N\ < N. 

To close the algebra one has to include all anti-commutators between supercharges, 
i. e. the full algebra based on two pairs of supercharges and P ± with real intertwining 
operators. Two supercharges generate two Polynomial SUSY, 

{K + ,K-} =V N (H), {P + ,P-}=V Ni (H). (121) 

The closure of the extended, Af = 2 SUSY algebra is given by 

{ K-,P+ } s ( 122 ) 

Evidently the components of operators 7Z, TV are differential operators of iV + N\ 
order commuting with the Hamiltonians h^°\ hence they form symmetry operators 
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K, TV for the super-Hamiltonian. However, in general, they are not polynomials of the 
Hamiltonians and these symmetries impose certain constraints on potentials. 

All four operators Vn(H), Vn^H), 71, TV commute each to other. The Hermitian 
matrix characterizing this M = 2 SUSY, 



is degenerate. Therefore it seems that the two supercharges are not independent and 
by their redefinition (unitary rotation) one might reduce the extended SUSY to an 
ordinary M = 1 one. However such rotations cannot be global (constant) and must use 
non-polynomial, pseudo-differential operators as "parameters". Indeed, the operator 
components of the "central charge" matrix 2(H) have different order in derivatives. 
Thus, globally the extended nonlinear SUSY deals with two sets of supercharges but 
when they act on a given eigenfunction of the super-Hamiltonian H one could, in 
principle, perform the energy-dependent rotation and eliminate a pair of supercharges. 
Nevertheless this reduction can be possible only after the constraints on potentials have 
been resolved. 

Let us find the formal relation between the symmetry operators 1Z, TV and the 
super-Hamiltonian. These operators can be decomposed into a hermitian and an anti- 
hermitian parts, 



The operator B is a differential operator of even order and therefore it may be a 
polynomial of the super-Hamiltonian H. But if the operator 8 does not vanish identically 
it is a differential operator of odd order and cannot be realized by a polynomial of H. 

It can be proven [69] that the hermitian operator B is a polynomial of the super- 
Hamiltonian of the order Nj, < N — 1. Let's use it to unravel the super-Hamiltonian 
content of the operator £, 



which follows directly from (11231) and (I124p . As the (nontrivial) operator S(H) is a 
differential operator of odd order N e it may have only a realization non-polynomial in 
H being a square root of (I125P in an operator sense. This operator is certainly non- 
trivial if the sum of orders N + Ni of the operators and is odd and therefore 



The existence of a nontrivial symmetry operator 8 commuting with the super- 
Hamiltonian results in common eigenstates which however are not necessarily physical 
wave functions. In general they may be combinations of two solutions of the Shrodinger 
equation with a given energy, the physical and unphysical ones. But if the symmetry 
operator £ is anti-Hermitian in the Hilbert space spanned on the eigenfunctions of 
the super-Hamiltonian H then both operators have a common set of physical wave 
functions. This fact imposes quite rigid conditions on partner potentials yW'W. It can 




det [2(H)] = V N V Nl - TLTV = 0,(123) 




S 2 (H) = V N (H)V Nl (H)-B 2 (H), 



(125) 



N e = N + N X . 
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be proven [69] that they inevitably belong to the class of transparent or reflectionless 
ones |123l 1124] . Such a symmetry has relations to the Lax method in the soliton theory 
[261127]. 

As the symmetry operator £ is anti-Hermitian its eigenvalues are imaginary but, 
by construction, its coefficients are real. Since the wave functions of bound states of the 
system H can be always chosen real we conclude that they must be zero-modes of the 
operator £(H), 

5(3)^ = ^ = 0, V N (E i )V Nl (E i )-B 2 (E l )=0, (126) 

which represents the algebraic equation on bound state energies of a system possessing 
two supersymmetries. Among solutions of (11261) one reveals also a zero-energy state at 
the bottom of continuum spectrum. 

6. Shape invariance for SUSY related potentials 

6.1. Shape invariance in linear SUSY 

One of the goal of SUSY QM design consists in search for exactly solvable models 
of Quantum Mechanics. For one- dimensional Schrodinger equation, a list of exactly 
solvable is well known: harmonic oscillator, Coulomb, Morse, Poschl- Teller, Scarf, 
Eckart potentials etc |126j . Each of such potentials has its own history of solution, 
but all of them were reproduced in the framework of a single algebraic procedure of 
Factorization Method |33j in the middle of last century. The related method was 
formulated in the context of modern Supersymmetric Quantum Mechanics with the 
help of new notion - shape invariance - introduced by L.E.Gendenshtein in [43j. 

The basic steps of standard shape invariance in absence of spontaneous SUSY 
breaking [13], |127] - |129j . jl] are the following. Let us consider a parametric family 
of one- dimensional superpartners h^(a), h^(a) and first order supercharges q ± (a) 
depending on a parameter a. We say that the Hamiltonians are shape invariant when 
in addition to supersymmetrical intertwining relations (j4j) they have the property 
(analogous to "cyclic" in [TTHIITT2 J). 

hF>(a) = h w (a)+Tl(a) } (127) 

where a = a(a) is a modified value of parameter depending on a and 71(a) is a (c- 
number) function of a taken positive for definiteness. Then the absence of spontaneous 
breaking of supersymmetry for all values of a implies that just the lowest eigenvalue 
Eq(o) of h(°\a) vanishes, and the corresponding eigenfunctions ^ (a) are normalizable 
zero modes of q~(a). This property allows to solve the models (i.e. to find the spectrum 
and all bound state wave functions) algebraically. 
To perform it we start from 

/i (0) (aK 0) (a) = E (a)^\a) = 0. (128) 
Consider the relation (j!27p to obtain, 

h {1 \a)^\a) = 7Z(a)^\a). (129) 
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We notice that ^/^(a) = ( a ) nas no nodes and therefore is the ground state wave 
function of h^(a). The combination of (J4]) and fl 1 2 9 [) yields, 

h^(a)q + (a)^\a) = TZ(a)q + (a)^\a). (130) 

Provided q + (a)^^\a) is normalizable, we have generated an excited state of h^°\a). 
These steps can be repeated up to the last step, where the resulting wave function \l> 
will not be normalizable anymore. There are notorious cases (oscillator-like potentials) 
where the spectrum is not bounded from above. The simplest and famous case is the 
harmonic oscillator with q = =f<9 + \x and lZ(a) = 2A. By the cyclic algorithm one 
easily reconstructs the energy spectrum En = 2n\. 

It is clear that the isospectrality of h^(a) and h^(a) (up to the only zero mode 
(o)) implies that there is no eigenvalue of h^°\a) between zero and the ground state 
energy EQ 1 \a) of hS l \ This observation can be used to prove that after a number of 
iterations one gets the entire spectrum of h^\a). Thus the shape invariance or cyclic 
method leads to spectrum generating algebras which allow to find algebraically all energy 
eigenvalues and eigenfunctions of shape invariant Hamiltonians in one-dimensional 
SUSY QM - complete solvability. 

We notice that the shape invariance relation (11271) combined with SUSY 
factorization of Hamiltonians ([3]) entails the conventional Heisenberg algebra, 

[q-,q + ]=K(a), (131) 

to supply q ± with the meaning of creation and annihilation operators. 

6.2. Intertwining with shift and higher order shape-invariance 

The particular construction of higher order shape-invariance can be realized [73] by 
breaking of hidden symmetry (11201) for two Polynomial SUSY algebras ( 11211) , 

h^k+=k+h^; k + N k- N = V N (h^)- k^k+=V N {hV>)i (132) 
h {0) P + Nl = P + Nl (h (1) + 2A); p + Ni p Ni = V Nl {h { % p~ Nl p + Nl = V Nl (h^ + 2A); (133) 

with the help of shift on a positive constant TZ(a) = 2A. In the sector of one of partner 
Hamiltonians (for example, h^) we obtain, 

/i(°) a + = a+(/i (0) + 2A), a-/i (0) = (/i (0) + 2A)a~ (134) 

if introducing the product operators a + = p^k^ = ( a )^ ■ We will call such a 
Hamiltonian as "higher order shape-invariant" . One can also work the other way around 
and start from ( 11 34ft with a differential operator a + of N\ + iV order in derivatives, 
factorize the latter as a product kjj with nonsingular operator multipliers and find 
an auxiliary super-partner Hamiltonian to obtain f l 1 3 2 D . 

Eq. (11341) is a spectrum generating, ladder ("dressing chain" |1101 1112j ) equation 
where a + plays the role of generalized creation operator which provides an excitation 
energy of 2A. In order to study the spectrum it is essential to determine normalizable 
zero modes of a - and a + . The former ones describe the lowest lying levels of the system 
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for positive 2A and one has to apply recursively the operator a + to them in order to 
generate the excitation spectrum. The energies of the zero modes can be obtained from 
vanishing the average of the operator product a + a~. In Nonlinear (Polynomial) SUSY 
this product can be evaluated algebraically because [45] 

a+a = p^k+p^ = P tV N (h^)p Ni = p N (h® - 2\)V Nl (hW), (135) 

see (11211) . In contrast to the simple harmonic oscillator, one has also the possibility of 
zero modes of the operator a + . Accordingly in this case the relevant operator product 
reads 

a -a+ = V N {h^)V Nl {h^ + 2A), (136) 

where the polynomials contain zeroes some of which correspond to zero modes of the 
operator a + . Thus ( 1134)) provides a connection between different levels of by a given 
shift, which is the simplest realization of the notion of shape-invariance [43]. Though 
some properties of the spectrum, like normalizable zero modes of a ± , will depend on the 
explicit product structure of a 1 * 1 , the equidistant excitation spectrum and corresponding 
wave functions can be mainly obtained algebraically . 

The generalized (or deformed) Heisenberg algebra is finally built by the following 
closure, 

[a+ a-] = V N (h® - 2\)V Nl {h {Q) ) - V N (h^)V Nl (h^ + 2A) ~ 2XV N+Nl - 1 (h^). (137) 

It is possible to study the certain consequences of ( I134p algebraically, using the deformed 
Heisenberg algebra ( 11371) without taking into account the specific definition of the 
operators a ± although the physical spectrum is essentially based on normalizable zero 
modes of a ± . 

The important comment is in due. Taking into account the minimization recipe 
from Subsection 3.4 and the Theorem 2 from Subsection 4.1 one can find in certain cases 
that the order of operators p Ni and k N can be lowered without changing the potential in 
the Hamiltonian h^°\ However the appropriate factorization of redundant polynomials 
of the Hamiltonian can be typically realized only at the level of individual SUSY algebras 
and, in general, cannot be extended to the deformed Heisenberg algebra ( I137j) . Thus at 
the formal level there may be several algebras ( I137P with the same spectrum pattern. 

6.3. Intertwining with shift and second order shape-invariance 

Let's elucidate the particular case when a 1 * 1 is of second order, in which we can obtain 
a generalized singular harmonic oscillator which is also shape- invariant. Correspondingly 
we consider in ( 11321) . ( 11331) k ,p of first order. We can thereby explore that they 
can lead to nontrivial consequences. The superpotentials which solve (1 132j) . ( II 33[) are 
growing linearly for large x and behave like 1/x at the origin just leading to two singular 
potentials pi] V i0 \x) and U (1) (x) as follows, 

v"> W = ^ + *£ + H, + i); 0«M = ^ + XV + MP - i)-(138) 
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These potentials are shape- invariant in the standard sense [43] and belong to the class 
of algebraically solvable models, because V^'(x; p, A) = V^°'(x; p + 1, A) — 2A. On the 
half-line with suitable boundary conditions for wave functions [45] at the origin the 
above potentials can be interpreted as radial harmonic oscillators for integer values of 
p. However this system makes sense also for arbitrary real p. 

The algebraic properties of these systems are based on the following factorization, 

h {0) = k + k~ =p + p- + \(2p + l); h {1) = k-k + =p-p + + \(2p-l), (139) 

satisfying the relations f)132p . fl 133[) . The intertwining operators are given by, 

k + = (k-y = -d-^-^; p + = -d-^ + ^. (140) 

In terms of the product operators a + = p + k~ and a = k + p~ one obtains an algebra 
suggestive of a generalization of the standard harmonic oscillator algebra, 

[h®,**] = ±2Aa ± , (141) 

where the raising/lowering ("creation/annihilation") operators a 1 * 1 can be presented in 
the form convenient for further evaluations, 

a+ = (a"). = exp {^} (-* + ^il) exp {-^} . (142) 

Similar results hold for hP-\ 

Meanwhile the crucial polynomial algebra describes the individual products, 

a+a = (h {0) -2A)(/i (0) - A(2p + 1)); a a+ = /i (0) (/i (0) - A(2p - 1)). (143) 

The zeroes of polynomials in f!143j) may indicate the zero modes ("vacuum states") of 
the annihilation, a~ or creation, a + operators. In the case under discussion it happens 
that the appropriate normalizable zero modes appear for the operators k + and p~ , 

f Ax 2 ] _ ~ ~ f Xx 2 ) 

& V>o,i = 0, ip ,i =x p exp < — — > ; p Vo,2 = 0, Vo,2 = x p exp < — — > , (144) 

which entails the two zero modes for the annihilation operator, 

f Xx 2 ) { Xx 2 ) 
a il) 0ti = 0; ^o,i = x p exp I — \ ; ^0,2 = x 1 p exp I — \ . (145) 

Accordingly they are associated to two zeroes of the first polynomial in (I143p . 

(h^ - 2A)Vo,i = = (/i (0) - A(2p + 1))^ 0>2 . (146) 

Another couple of zeroes for the second polynomial in ( I143j) are not related to any 
normalizable solution. 

Eventually the oscillator algebra f!14ip generates one ladder of equidistant levels 
for p > | or, equivalently for p < — | and two independent ladders of equidistant levels 
for — | < p < | or, equivalently, ~ < p < | . For p = | these two ladders coincide. 
The corresponding spectra of are: E n i = 2X{n + 1) or E n>2 = X(2p + 1 + 2n) for 
n = 0,l,.... 
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The generalized Heisenberg algebra for this sort of oscillator reads, 

[a+a-] = -4A/i (0) + 2A 2 (2p + l). (147) 

There might be an (superficial) impression of the algebra f l 1 4 1 j) as to be a broken 
hidden symmetry of the M = 2 SUSY (I12ip . (j!22p . However the contraction A — > 
results in the system with a ± ~ and continuum spectrum and therefore no any 
nontrivial hidden symmetry emerges. 

A nontrivial, shape-invariant, breaking of the M = 2 SUSY arises for the third-order 
spectrum-generating algebra. 

6.4- Intertwining with shift and third order shape-invariance 

Let us now consider the intertwining operator to be a reducible or irreducible 
operator of second order but p ± still of first order like in (I104H [73], |130j . Accordingly 
they are parameterized by 

^ = d 2 T2f(x)d + b(x)T f(x); (148) 

p ± = T d + x{x). (149) 

The solutions of (11321) are unchanged in respect to Eqs. (11051) of Subsection 5.1 but 
Eq. (ll33p implies a shift for the potential V^(x) in (11061) . Thus the consistency equations 
for xi x ) an d f{x) are modified. The first Eq. (11081) becomes: 

X (x) = 2f(x) + Ax, (150) 

where a possible integration constant can be ignored because of a shift of x, which fixes 
the origin of coordinate. Then the potential V^(x) can be written from (11051) as: 

V {0) {x) = -2f(x) + Af 2 {x) + AXxf{x) + X 2 x 2 - X. (151) 

The coefficient in the intertwining operator k reads, 

b = - (2f 2 + 4Xxf + 4AV + a) . (152) 

In the second Eq. ( 11081) one has to replace \ 2 in accordance with ( 11501) . i.e. f[x) satisfies 
the following equation, 

/" = Qr\ + 6 / 3 (*) + 8Aa;/ 2 (a;) + 2(AV - A + a)f(x) + (153) 

The equation (11531) can be transformed by the substitution f(x) = l/2y/Xg(y); y = \f~Xx 
to the Painleve-IV equation [74] , 

9 " = CM + ^ 3{y) + Ayg2{y) + 2{y2 " a)9{y) ~ Wi' (154) 

where 

a , 4c? , 

a = l-^; b = -—. (155) 

This equation has been studied intensively in the last years [131j . [T3~2"] .In what follows 
we will focus mainly on asymptotic properties of its solutions which will determine the 
asymptotics of potentials ( I15ip and the normalizability of eigenfunctions. 
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The spectrum generating operators a + = p + k = (a y form the polynomial 
algebra, 

a+a = /i (0) [(/i (0) + a - 2A) 2 + d] ; a a+ = [/i (0) +2A] [(/i {0) + a) 2 + d] .(156) 

We show how one can derive the spectrum from fjl34j) and (11591) if normalizable zero 
modes of the annihilation operator a - exist. We stress that this algebraic method is 
very powerful now since the explicit form of the potential is known only in terms of 
Painleve transcendents. 

The equation for zero modes of a - reads: 

a -^[ 0) = k+p'^f = (157) 

where k labels the normalizable solutions. 

In accordance to the non-linear SUSY algebra (I156p . the equation for eigenvalues 

4 0) is, 

£(0) . [(£(0) + a _ 2A) 2 + d\ = 0, (158) 

and has at most three real solutions. 

Concerning the generalized Heisenberg algebra |133] , the modification in respect to 
(I14ip in the previous Subsection is given by, 

[a+, a"] = -2A(3(/i (0) ) 2 + (4a - 2A)/i (0) + a 2 + d) (159) 

and similarly for h^ 1 '. 

6.5. Spectrum patterns for third order shape-invariance 

We give now a short description of typical spectrum patterns (the details can be found 
in [T3]), [130J. Below on, the case of equal asymptotics of f(x) at ±oo is discussed. 

• Three normalizable zero modes of a~ and, respectively, three equidistant sequences 
of levels may arise only for A > \/—d if —y—d > a > y 1 — d — 2 A or 
2A — \J—d > a > y—d. The corresponding solutions of Painleve-IV equation 
must have the leading asymptotics f(x) ~ —Ax/3 with sub leading oscillations 

cos(Xx 2 /V3 + o(x)). 

• Two normalizable zero modes of a - (and two sequences of levels) may exist for 
different solutions of Painleve-IV equation. 

— Namely, for A > y—d the solution with asymptotics f(x) ~ —Ax/3 provides 
the fall off of * W (x), $f(x) if a < y^d - 2 A and of ^{x), * { °\x) if 
2A + y^d > a > 2A - y^d. 

— For A > y 1 — d the solution with asymptotics f(x) ~ —y / —d/2Xx generates 
two sequences of levels starting from \&q (x), ^_(x) if 2Xn + y 1 — d < —a < 
2A(n + l) - V^d; n = 0,1,2... 

— For A < y 1 — d two sequences (from ^/^(x) and \l/^(x)) are generated for two 
possible asymptotics f(x) ~ y/^d/2Xx and f(x) ~ —Ax/3 if 2A + y'—d > a > 
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For a > 2X + y/—d and arbitrary positive A one sequence of equidistant levels will be 
realized with ground state \&q (x) and one of three asymptotics: f(x) ~ —As/3, ~ 
±v /z c!/2Ax. 

— For f(x) ~ — Ax/3 only one sequence of levels will be realized also for 
A = y/—d, a = ±y/—d (starting from (a;)), for A = y/—d, a < \J — d 
(starting from ^^\x)) and for A < y/—d, a < 2A + y/—d (starting from 
*f(*)). 

— For asymptotic behavior f(x) ~ — \ / —d/2Xx and A > a 
consists of one sequence of levels if \a + 2Xn\ < V~d, n 
from ty_\x)) or if a = — \/— d (starts from ^^)- 

— Last, for f(x) ~ \/—d/2Xx one obtains one sequence with ground state \P+ if 
A < \f--d and a < y/—d or if A > y/— and |a + 2An| < y/—d, n = 0, 1, 2, ... 

For specific values of parameter a the normalizable zero mode of a + can cause 



-d the spectrum 
0, 1, 2, ... (starts 



(°), 



x) 



<0) 



0,£ {0) <0, 



truncation of one of two sequences ( with ground state 
Let us also remark that particular value a = — 2A — y/—d, i.e. E\ 
just corresponds to the case when Painleve-IV equation (I154p has a class of 
particular solutions which coincide with solutions of the Riccati equation g'(y) = 
9 2 (y) + 2y9(y) + Vb- Substituting this Riccati equation into Eq. (I151j) one finds that 
potential is the pure harmonic oscillator: V\{x) = X 2 x 2 — A. 
One additional singlet state satisfies the equation 
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The doublet representation (\I/q (x) 
f)134p . (1159P is built on solutions of the equation: 



—d/2Xx. The spacing between the two ground 
2A =p 2v^d. 

3l + ^^\x)) of the spectrum generating algebra 
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a-¥°\x) 



0. 



It may hold when a 



-2A 



-d for arbitrary positive value of A and when 



a = —2X + v ' ~d for A > \J—d. It is equivalent to a+^Q^x) = \l/ li2 (x), which is 
satisfied when f(x) obeys the following equation: 

8A/ 2 (x)(f (x) + 2/ 2 (x) + 2Xxf(x) — 2A =F V^d) = (/'(x) + 2f 2 (x) + 

2Xxf(x) t \^d) ■ \ (f'(x) + 2f 2 (x) + 2Xxf(x) - 2A T v 73 ^) 

(-/'(x) + 2/ 2 (x) + 2Ax/(x) — 2A =f v^) - 4A(A ± V^d) 

One can show that all solutions of this equation fulfill the Painleve-IV equation 
f)153p . These solutions have the asymptotics f(x) ~ ±\ / —d/2Xx and cannot have 
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any (pole) singularity. The spectrum consists of a doublet (0, 2A) and infinite 
sequence E n = ±2 V /Z d + 2(n + 2)A, n = 0, 1, 2... 

Recently some higher-order generalizations of shape-invariance have been elaborated 
leading to the Painleve-V etc. equations [134] . [T35] . As well new solutions of Painleve- 
IV with complex parameters based on the third-order shape invariance have been found 



7. Nonlinear SUSY for non-stationary Schrodinger equations 

7.1. Linear SUSY and hidden symmetry 

In this Section our aim is to elucidate that many of the nonlinear SUSY constructions 
illustrated before can be implemented also in the Schrodinger time dependent framework 
[137], [75], [6], [138]. 

The first- and higher-order intertwining operators and the corresponding relations 
between non-stationary one-dimensional Schrodinger operators can be introduced 
straightforwardly. But, as compared to the stationary case, already the first-order 
intertwining relations imply some hidden symmetry [§| which leads to a specific quantum 
dynamics when the evolution is described by quantum orbits and results in the R- 
separation of variables [75J. In turn second-order intertwining operators |75]. |139j and 
the corresponding non-linear SUSY give rise to the quantum motion governed by the 
spectrum generating algebras. We start with the non- stationary Schrodinger operator 

S[V] =id t + dl-V(x,t) . (162) 

Here d t = d/dt and d x = d/dx denote the partial derivatives with respect to time and 
position. When they are applied to some function /, the following notations for these 
derivatives are used: f(x,t) = (d t f)(x,t) and f'(x,t) = (d x f)(x,t). 
The intertwining operator of first order [75] is given by 

q+ = ^ {x,t)d t + ^{x,t)d x +^{x,t) (163) 

with, in general, complex-valued functions £o>£i an d ^2-The possibility of a 
complexification of the intertwining (Darboux) (first and also higher order) operator 
was emphasized by [139J. 

For the Schrodinger operator (I162p the intertwining relation reads 

S\V<®]qt = qtS[V^] , (164) 

where the functions & (i = 0, 1, 2) and are not independent of each other. It can 

be also represented in the SUSY form Eqs. ( jHl) .(fl2l) when the stationary Hamiltonians 
are extended to the Schrodinger operators iSfV^ ^ 1 )], then 

§ Symmetry properties of the time dependent Schrodinger equation were studied in general in [140) . 
[141) (and references therein). 
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After inserting (I162p and the intertwining operator (II 6 3 p into relation (11641) it can 
be derived that £ and £i may depend only on time, i.e. £q = = £[. Then the 
assumption that £o does not vanish identically would entail that the potential difference 
y(o) _ yW does depend on time only. This is a rather trivial case and therefore let 
us take £o = . Making the appropriate choice of variables = e % P® p(t) and 

£,2(%,t) — e 1 ^^ p(t)tu'(x, t) with real (3, positive p and complex to functions one finds 

V(°) (x, t) = u' 2 (x, t) + lu"(x, t) - uj(x, t) + ait) - Pit) + ipit) I pit) , 

(166) 

V w (x, t) = tu' 2 (x, t) - u"(x, t) - ico(x, t) + a(t) , 

where a is a time-dependent complex- valued integration constant. One may set [73] 
(3 = without loss of generality. Furthermore, one may also remove a — > as being 
absorbed in u by the shift u> — > u> — i J dta. Finally, we are left with 

„( I , 1) = l /.( I ,«)^( I ,«)-M, 1 ) + «.)M.), 

(x, t) = co' 2 (x, t) - u"(x, t) - ico(x, t) . 

Here the super-potential u is not arbitrary as the potentials are assumed to be real. It 
could be achieved by taking a stationary real super-potential. However this choice leads 
to the stationary SUSY QM discussed previously . Therefore, we consider a complex 
super-potential 

cj(x, t) = u R (x, t) + icui(x, t) (168) 

with real functions co R and ujj. The reality condition Im = Im = is fulfilled 
if 

2(tu R )" + p/p = , 2{u R )'{u I )' - {ui)" -co R = 0, (169) 
which can be integrated to 

wiOM) = -]^x 2 + \ P (t)p(t)x + 1 (t) , 

4 PW 2 (170) 
u R {x,t) = - lnp(t) + K(x/p{t) + //(*)) , 

where p and 7 are arbitrary real functions of time and K is an arbitrary real function of 
the variable y = xj p + p. In terms of these functions the final form of the two partner 
potentials is 

VV>»Hx,t) = -L- ) [K a (y)±K>>(y)] 

m x > + (mm + p M) * - pHmt) + m m 



4p(t) V 2 / 4 

and the intertwining operator reads 



q+(x, t) = p(t)d x + K'(x/p(t) + /x(t)) - % - (p(t)x - p 2 (t)p(t)) ■ (172) 
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Let us demonstrate that the non- stationary Schrodinger equation S[V^°' ) '^}ip^ oS> '^ = 
with potentials given in Eq. (I17ip (which is equivalent to the intertwining (I164p ) admits 
a separation of variables. Indeed, after the transformation 

y = x/pQ+tit), Tp mi) (x,t) = -^ e -^^(f> mi) (y,t) = n(x,t)(f>^(y,t)(m) 

this Schrodinger equation becomes quasi-stationary [142J 

ip 2 (t)d t <j>^\y,t) = [-d 2 + K'\y) ± K"{y)] ^(y,t) , (174) 

which is obviously separable in y and t. Hence, the solutions of the original 
Schrodinger equations have the general form i()(x,t) = Q(y , t)Y (y)T (t) which is known 
as the .R-separation of variables [143] . In other words, for any pair of Schrodinger 
operators 5[V^^ - ) ' < - 1 ^, which admits a first-order intertwining relation (I164p there exists a 
transformation (11731) to some new coordinate in which the potentials become stationary 
(see also |138j ). 

This /^-separation of variables is certainly related to the existence of a symmetry 
operator for the super-Hamiltonian. First, one can directly verify the adjoint 
intertwining relation for real potentials, 

<kS\y<®) = s\y<»]qr (175) 

where 

«T = = ~P(t)d* + K'(x/p(t) + /i(t)) + i (p(t)x - P 2 (t)p(t)) ■ (176) 
Then from (11641) . (11651) and (j!75j) we obtain the closure of the SUSY algebra, 

{Qt,Qt}=n t , [4t^]=o, Qt = (QtV, (177) 

where the symmetry operator 1Z t has the following components 
Rf = qfqf = -p\t)dl + 1 {pp(t){x,d x } - 2p 3 (t)p(t)d x ) 



+ 



K'(x/p{t) + p(t))] 2 ± K"(x/p{t) + p{t)) + - A {p{t)x - p\mt))- 



exp{-iwi(s,*)} [-d 2 y + K' 2 {y) ± K"{y)] exp{iu T (x,t)}. (178) 

Thus the quasi-stationary Hamiltonians in Eq. (11741) are just unitary equivalent to the 
symmetry operators Rf. It means that the supersymmetry entails the separation 
of variables because it provides a new symmetry for the super-Hamiltonian. As 
consequence the quantum dynamics splits in orbits with a given eigenvalue of the 
symmetry operator. 

7.2. Second- order intertwining for stationary potentials: spectrum generating algebra 

Let us relate a pair of Schrodinger operators S[Vi] and <S[14] by second-order 
(intertwining) operators of the form, 

qf(x, t) = G(x, t)d 2 x - 2F(x, t)d x + B(x, t) . (179) 
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We will explore the existence of the time-dependent SUSY charges with appearance of 
the spectrum generating (oscillator like) algebras for corresponding Hamiltonians. 

As in the first-order case it can be shown [75] that the inclusion of an additional term 
of first order in dt leads to the trivial solutions when the difference V\ — V% depends on 
the time t only. Furthermore, from the intertwining relation (11641) with the above one 
can conclude that the function G cannot depend on x and it is even possible to exclude 
a phase. In other words, without loss of generality G(x,t) = g{t) and accordingly, from 
now on, an intertwining operator can be reduced to, 

q+(x, t) = g(t)d 2 x - 2F{x, t)d x + B(x, t) . (180) 

In [75] particular solutions of the intertwining relation ( 11641) were constructed with qf 
as given above . In this section we shall analyze the solutions of the intertwining relation 
(11641) for the case where both potentials V\ and V 2 are stationary, i.e. do not depend 
on t. 

One class of such solutions is known from [45J. Assuming a supercharge qf with 
real coefficient functions independent on t, one finds that the corresponding solutions of 
(11641) coincide with those of the stationary intertwining relations (— d 2 + V^°\x))q + {x) = 
q + {x)(-d 2 x + V^{x)). 

Here we are interested in more general solutions of ( 1 164ft when operators qf depend 

on t, 

{id t - h^)q+{x,t) = qt(x,t)(id t - h,W) , (181) 

with stationary Hamiltonians ft,(°)>W = — d 2 +V^'^(x) but time-dependent intertwining 
operators. 

Let us employ the suitable ansatz with simple t-dependence in (11801) . 

q+(x, t) = k + (x) + A(t)p + (x) , (182) 

where, 

k + (x) =d 2 x - 2f(x)d x + ~b(x) - f'(x), p + (x) = -d x + X (x). (183) 

Here all functions besides A are considered to be real. We also assume A ^ 0. With 
this ansatz the intertwining relation (11811) can be shown [75j to yield, 

iA = 2rh + 2mA , 

k+h^ - h^k + - = 2 mp + , (184) 

h(P) p + -p+hM = 2mp + , 

with real constants m and m. 

We find it interesting to focus on the case m ^ to explore certain spectrum 
generating algebras. The first equation in ( 1184jl immediately leads to 

A(t) = m e' 2imt - rh/m (185) 

with a real mo and, 

/ 771 \ ~ 777 

q+{x, t)=d 2 x - (2f(x) + -)d x + b(x) - f'{x) + - X (x) - m e- 2mt p + (x) . (186) 

\ 777 / 777 
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Without loss of generality we may set rh = because a non-vanishing rh may always be 
absorbed via a proper redefinition of / and b, i.e. of the operator k + . 

As a consequence, the second relation in fjl84j) leads to a second-order intertwining 
between and ft/ 1 -*. This has already been considered in [45J and it was found that 
the potentials V(°\ and the function b can be expressed in terms of / and two 
arbitrary real constants a and d as in Subsect.5.1, Eq. fll05[) . 

f"{x) f\x) d 



V mi) (x) = T 2f(x) + f(x) + 



2m 4P(x) AP{x) ' {m 



b(x) = -f{x) + f(x) - J —^r + + 



2f(x) AP{x) 4/2(x) ■ 
The corresponding second-order SUSY algebra generated by the supercharge K is similar 
to ffTOTl) . 

One may find some similarities between the present intertwining algebra ( I184p and 
the extended SUSY relations discussed in Section 5. But we emphasize that now for 
m 7^ the last relation in ( 1184p does not generate a second SUSY. Rather it creates 
the equivalence of relatively shifted spectra of two Hamiltonians h + and h~ which is 
typical for spectrum generating algebras emerging for the shape-invariant Hamiltonians 
(Sect. 6). Specifically, 

p + p~ = — m + c; p~p + = hP~> + m + c, (188) 

where c is a real constant. Therefore the reflectionless potentials found in Section 5 are 
produced only in the limit of m — 0. 

The genuine spectrum generating algebra for stationary Hamiltonians h^ '^ 1 ' can 
be derived from Eq. (1 184ft 

l hi ° ] i a (o)l = ~ 2ma (o) ' >) = k+p ~ ' 
a+ ] = 2ma+ 0) , a+ 0) = p+k~ , 

a (i)] = 2ma (i) > a (i) = ' 

a 5)l = ~ 2ma t) > a a) = p~ k+ > 

where p~ = (p + )^, = (& + )t and a^ ^ = (a^^)*. 

The closure of this spectrum generating algebra is a polynomial deformation of 
Heisenberg algebra [133] (see Subsect. 6.4), 

HoAV*m}= pm (h^). (190) 

The explicit form of the polynomials p( 0,1 can be obtained with the help of Eqs. (1 184ft 
and (1188ft for rh = 0. For instance, 



a; nN a 



(0)~(0) 



— m + c) \(h i0) - 2m + a) 2 + d 



a- 0) a+ 0) = (h® +m + c) (/i (0) + a) 2 + d , (191) 
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where the notations from Eqs. (11071) and (j!88p are employed. The polynomials P^ ' 1 ) 
turn out to be different for the isospectral partners /i ± , 

p(°\h {0) ) = -6m(/i (0) ) 2 + 4m(2m - 2a - c)/i (0) - 2m[a 2 + d + 2(a-m)(c- m)]; 

= _6 m (/ l ( 1 )) 2 - 4m(2m + 2a + c)/i {1) - 2m[a 2 + d + 2 (a + m)(c + m)]. (192) 

Hence the two spectrum generating algebras are, in general, different that is essentially 
due to the shift in intertwining relations (1188!) . When comparing to the equations in 
Subsection 6.4 we find that the generating algebras (I19ip . (11921) are in fact originating 
from the third-order shape invariance. 

There is a formal discrete symmetry between their constants and Hamiltonians 
hW,a,c=>-hW,-a,-c. 

The intertwining relation (I18ip and its adjoint give rise to the symmetry operators 
qtqt and g t ~g t + for (id t -h^) and (id t -h^), respectively. Using Eqs. ffT07]) . ffT85|) .f rT88|) 
and after elimination of polynomials of the Hamiltonians h^ ' 1 ^ these symmetry operators 
may be reduced to, 

(193) 

As our potentials do not depend on time the time derivatives R^ 0,1 \x, t) of Hermitian 
symmetry operators i?^ 0,1 ^(x,t) form the independent set of hermitian symmetry 
operators which do not commute between themselves. Similar results have also been 
obtained in |141] using a different approach. We see that the non-stationary SUSY 
delivers the time-dependent hidden symmetry operators [14J, 

e2m %i) = 7T^ (0,1) (M) - -^—R^(x,t), (194) 
v ' ; 2mo Ammo 

which give the entire set of spectrum generating algebras previously found for the third- 
order shape invariance. Thus hidden symmetries for nonstationary Schrodinger super- 
Hamiltonians give the universal framework for both stationary hidden symmetries and 
for spectrum generating algebras of shape invariant SUSY systems. 

8. Polynomial SUSY QM in d = 2 

8.1. Supercharges of second order in derivatives: generalities 

It was shown in Section 2, that multidimensional generalization of SUSY QM includes 
both scalar and matrix Hamiltonians, and some physical examples of matrix problems, 
incorporated into SUSY QM, were given in Subsections 2.2 and 2.5. In general, 
the construction does not provide any intertwining and therefore any direct relations 
between spectra and wave functions of two scalar components of super-Hamiltonian, if 
we are not interested in models with separation of variables. But such opportunity can 
be opened by applying and development of the ideas of polynomial SUSY of Section 3 for 
higher dimensionality of space. In particular, this approach turned out [76| - [80], |82j, 



R (1 \x,t) = m 



„2imtr.— _i_ „— 2imt c ,+ 

e a (1) -+- e a (1) 
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[83] . [85] - [88], [90] - [93J to be very fruitful for study of nontrivial spectral problems 
in the case of d — 2 (see also [T3], [16J). 

Thus, in order to get rid of matrix component of two-dimensional super- 
Hamiltonian, we shall explore the supercharges of second order in partial derivatives. 
The simplest variant of second order supercharges - of reducible form q ± = {qf){qf) - 
is not very promising: the intertwined partner Hamiltonians differ by a trivial constant 
only, and both of them admit the separation of variables (see details in [TB], [77]). By 
this reason, here we are interested in general irreducible second order components of 
supercharges, 

q + = g lk {x)did k + Ci{x)di + B(x); q- = (g+) f , (195) 

with real "metrics" gi k (x) and coefficient functions Ci(x), B(x). The familiar intertwining 
relations 

h (o) q + = q + h (i). q - h (o) = h (i) q -_ ( 19 g) 
for two scalar two-dimensional Hamiltonians of Schrodinger form, 

/l (0),(l) = _ A (2) + v (0),(l) (f) (1Q7) 

can be rewritten as a system of nonlinear partial differential equations for gi k (x), the 
coefficient functions Ci(x), B{x) and potentials V^°'^(x). The general solution for the 
"metrics" gi k (x) is [76], [77], 

gn = axl + a^ 2 + h; #22 = ax\ + a 2 x x + b 2 ; #12 = -\\2ax\x 2 + a x xi + a 2 x 2 ) + b 3 , (198) 

where (a, a i; 6j = Const). Thus, the senior in derivative part of supercharges belongs to 
the E(2) - universal enveloping algebra [143] . We distinguish four different (inequivalent) 



classes in second derivatives [76], [77] . 

q W+ = 7 (P* + P*)+C i d i + B; (199) 

g( 2 )+ = a Pl + (a + 7 )Pf + CA + B; (200) 

q {3)+ = Pi} + a(P? + Pi) + CA + B; (201) 

= 7 J 3 2 + pP* + aPl + CA + B, (202) 



where J3 and Pi, P 2 are generators of rotations and translations, correspondingly, and 
7^0. 

The solutions 0198p must be inserted into six other equations which followed from 
the intertwining relations (I196p . 

diC k (x) + d k Ci{x) + A^g tk (x) - (V^(x) - VW(x))g ik (x) = 0; (203) 
A^C t (x) + 28 l B{x) + 2g lk {x)d k V ( - l \x) - (V^(x) - V^(x))C t (x) = 0; (204) 
A®B(x) + g ik {x)d k diV {l \x) + C l (x)d i V {1 \x) - (V^°\x) - V^(x))B(x) = 0. (205) 

It is clear that we have no chances to solve these equations with general form H198[) of 
gi k (x). Therefore, one is forced to consider the simplest cases of constant metrics g ik 
looking for particular solutions for functions Ci(x), B(x), V^°'(x), V^\x). 
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The first sample of constant metrics is the elliptic (Laplace) one corresponding to 
(I199p . g ik = 5 ik . In this case, the system (I203p - (I205p is essentially simplified, so that all 
coefficient functions can be found analytically [76], [77] . In particular, the combination 
C = C\ + 1C2 depends only on z = x\ + ix%, and it has the specific form, 

C 2 (z) = az 2 + bz + c, (206) 

where a is real and b, c complex constants. All other functions can be also found 
explicitly, but depending on values of constants in ( 12061) . both potentials allow 
R— separation |143] of variables - in elliptic (for b = 0, a 7^ 0), parabolic (for 
a = 0, b 7^ 0) or polar (t = c = 0, a / 0) coordinates. The situation with both 
a 7^ 0, b 7^ is reduced to the first variant by a suitable constant shift of coordinates. 
The case of polar coordinates just corresponds to reducible supercharges qfqf. Thus, 
for all possible values of coefficients the class of two-dimensional problems with elliptic 
metrics of supercharges can be reduced to two one-dimensional models, and it will not 
be considered below. 

8.2. Hyperbolic (Lorentz) metrics 

8.2.1. Supercharges with twisted reducibility. The constant metrics of hyperbolic 
(Lorentz) type, g ik = diag(+l, — 1) gives much more interesting models [77], [79], [75]. 
One specific class of such solutions of intertwining relations (11961) with the so called 
supercharges with twisted reducibility. 

q- = (g+)t = q+(a 3 ) lk q^ = + d t W {x)){a 3 ) lk {+d k + d k W(x)), (207) 

was studied in papers [82], [85], [87]. In (12071) . W(x),W(x) are two different functions 
(superpotentials), 03 is the Pauli matrix, and summation over k, I = 1, 2 is implied. This 
form of supercharges can be considered as a generalization of the simplest reducibility for 
elliptic metrics, mentioned in the beginning of previous Subsection. In contrast to that 
case, the form (I207j) leads to very nontrivial models, which do not allow the standard 
separation of variables. This construction of supercharges can be interpreted as obtained 
after gluing of two systems of two-dimensional SUSY QM of Subsection 2.1 with first 
order supercharges with different superpotentials W, W, correspondingly. These systems 
are glued by their matrix components of super-Hamiltonian: these components are taken 
coinciding but up to an unitary constant matrix rotation by 03. In such a case, both 
Hamiltonians = qfq^ and = q^q{ are quasifactorized in terms of qf and qf, 
correspondingly. 

It was shown in [82], [85] in a general form that Eqs. (jl96p and (1207P lead to the 
following representation for superpotentials W, W in terms of four functions /ii,2,±, 

X = (J.l(xi) + yU 2 (x 2 ) + fl+{x+) + 

X = /ii(si) +MZ2) - L l +( x +) ~ A*-0&-)> 
where the light cone variables x± — X\ ± xi were introduced. These functions /^i,2,± 
have to satisfy the following equation, 

Mi(zi) [n' + {x+) + + /i 2 (x 2 ) [fj! + (x+) - /x'_(a;_)] = 0, (208) 
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where fJ>'(x) means derivative over the argument. By substitutions (ft = /i', it is evident 
that we deal with a purely functional equation with no derivatives, 

Mm) [<!>+(*+) + <M*-)] = ~Mm) [Mm) - M*-)] • (209) 

Each solution of Eq. fl209|) gives the corresponding solution of the intertwining 
relations (I196P for the potentials V, V and for the supercharges q ± , 

v (0 \x) = (Mm) - Mm)) + (<t>l(x 2 ) - Mm)) + (f + {m) + V2<p' + ( x+ )) + 

+ (0 2 _(x_) + v / 20'_(x_)), (210) 

vU(x) = (Mm) - Mm)) + (Mm) - Mm)) + (f+(x+) - V2<p' + (m)) + 

+ (<P_(x-)-y/2<f/_(x-j), (211) 

q± = d\ - d 2 2 ± V2(Mm) + <j>-{x-))d x =f V2(^ + (x+) - <M^-))<9 2 - 
- (Mm) - Mm)) + (Mm) - Mm)) + 2Mm)M^~)- 

Thus, in order to find the systems with intertwining f)196p by supercharges of the form 
(12071) . it is necessary to solve (I209p . This equation seems to be rather complicated, but 
it appeared to be solvable in a general form (see details in [87J). In particular, it was 
shown that <p\^{x) are defined from solutions of the first order nonlinear differential 
equation, 

2 = a ^ 2 + bcf>l 2 + c, (212) 

where a, b, c are arbitrary real constants. All solutions of this equation can be expressed 
in terms of elliptic functions, and they are described for different ranges of parameters, 
for example, in Appendix B of [144] . In turn, the functions <j)± are defined uniquely 
by solution of (121 2ft : they also satisfy equation of the form (12121) . but with different 
constant coefficients. Depending on the values of a, b, c, solutions can be functions 
of their argument with finite or infinite period. Corresponding models were studied in 
detail in [87] (periodic potentials) and in [82], [85] (potentials with infinite period). 



8.2.2. Irreducible supercharges. Let us consider now the most general supercharges 
with hyperbolic (Lorentz) metrics without any kind of reducibility [76J, [77], [79], [13]. 
The light cone variables x± = x x ± x<i are again the most convenient to look for the 
solution of intertwining relations. In particular, it follows from the system ( 12031) - ( 12051) 
that the combinations C± = C\ =F C<i depend only on one argument: 

C+ = C + (x+), C_ = C_(x_). (213) 

Then, the whole system (12031) - (I205f) can be rewritten in a very compact form of two 
equations for functions C±(x±) and an auxiliary function F(x), 

04C-F) = -d + (C + F); (214) 
d 2 + F = dlF. (215) 
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It is easy to write the general solution F(x) of Eq. (1215j) : it is expressed in terms of two 
arbitrary functions, 

F = Fjx+ + xJ + F 2 (x + -xJ, (216) 

both are defined up to an arbitrary real constant, F\ — > F\ + 7, F2 — > F2 — 7. The 
potentials V^'^^x) and the function B(x) are expressed in terms of four functions 
Fj2xi), F 2 {2x 2 ) and C±(x±), which must satisfy the only equation (12141) . 

v mi) = ± l ( ^ + C ) + l (c a + ) + i (F 2 (x + - xJ) - F 1 (x + + *_)) ; (217) 
B = \{C+C- + F 1 (x + + xJ) + F 2 {x + - x-j), 

where C means derivative in its argument. Although Eq. (l214p has a compact form, it is 
very nontrivial being the functional-differential equation: the functions each depend on 
specific argument. As a rule, equations of this kind have no regular recipe for solution, 
and this is the case for Eq. (j214p . Usually, one may try to solve such equation starting 
from different suitable anzatses for functions C±(x±), F(x). 

1) Let us choose C_ = 0, then from (12141) one obtains F = <p(x^)/C + (x + ). After 
inserting into Eq. (l215j) the separation of variables is possible, and particular solution 
reads, 

C + (x + ) ' 



5i exp(\/A • x + ) + 5 2 exp(— y/\ ■ x + ) ' 
Fi j2 (2x) = 5icri t 2 exp(2VAx) + S 2 a 2: i exp(— 2VAx), 

where the Greek letters - arbitrary constants: depending on sign of A they may be 
real/ complex. 

2) Let us consider the anzats with factorized function F(x), F = F + (x + ) ■ 
Then from Eq.(EHD, 

x± 

c± = lt ± lkl F±dx>± ' (218) 

and there are two options to fulfill the condition (I215p . i.e. F(x) = Fi(2xx) + F 2 (2x 2 ), 

a) F±(x±) — e±x±, b) F± = a± exp(V / A ■ x±) + 5± exp(— vA ■ x±). (219) 

Corresponding potentials can be found according to Eq. (l217p . being similar to ones 
obtained in [I3_5J in quite different approach. 

3) Let us start now from the general solution of (12141) . 

F = L(! d 4±- f^=)/(C + C.). (220) 



C + J c 

Then Eq. (12151) gives the functional-differential equation for the functional L(A + — AJ) 
with A' ± = 1/C±(x±), 

AW Alll 

- -jf) L(A + - AJ) +3(Al + A"_)L'{A + - AJ) + {A% - A'*)L"(A + - A J = 0, (221) 
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where L' denotes the derivative of L with respect to its argument. If we take functions 
A± such that A"± = \ 2 A±, A = const, Eq.f l221j) will become ordinary differential equation 
for L with independent variable (A + — AS). It can be easily solved, 

L(A + - A_) = a(A + - A_y 2 + /3, 

where A± = a±exp(\x±) + 6±exp(—Xx±) with a + ■ 5 + = a_ ■ o~_ and a, (3 - real constants. 
For A 2 > 0, choosing a± = —6± = k/2 or a± = +5± = k/2, we obtain (up to an arbitrary 
shift in x±) two particular solutions, 

3a) A± = /csinh(Ax±), 36) A± = k cosh(Aa;±). 

Then (1220]) leads to, 
hi 

3a) F x (2x) = 27 — 7 + fc 2 cosh(2Ax); (222) 

cosh (Xx) 

F 2 (2x) = , - + k 2 cosh(2Ax); C± = . f. fc ^ 0, 

smh (Ax) cosh(Ax ± ) 



36) iq(2x) = -F 2 (2x) = . i 2 - + fc 2 sinh 2 (Ax), C± = fc ^ 0. (223) 



■ ,2/> \ + ^sinh 2 (Ax), C± = . f 
sinn (Ax) smh(Ax±) 

For A 2 < hyperbolic functions must be substituted by trigonometric ones. 

We have to remark that the case A 2 = 0, i.e. A± = 0, is not of interest, leading to 
trivial superpartners. However, choosing in H223I) A — > 0, k, ki, k 2 x — > simultaneously, 
so that A 2 ~ k\ ~ k 2 x ~ fc 2 , we obtain the solution, 

F 1 (2x) = -F 2 (2x) = k lX - 2 + k 2 x 2 , C± = — . (224) 

x± 

One can check that ( 1214ft is also satisfied by 

Fi(2x) = -F 2 (2x) = fcix 2 + fc 2 x 4 , C± = ± — . (225) 

x± 

4) Starting again from ( 12201) . it is convenient to pass on to new variable functions 
C± = ±/+//±- Then F in (12201) is represented in the form F = U(f + f-)f' + f'_ with an 
arbitrary!}]] function U. After substitution in ( 1215P one obtains the functional-differential 
equation, 

/ f " f"\ / f "' f"'\ 

U+f- - flf'- 2 )u"(f) + 3/(^± - ^) c/'(/) + ( - ^r) W) = o- / = /+/- 

For particular form of functions f± = a±exp(\x±) + f3±exp(—\x±), this equation 
becomes an ordinary differential equation for U with independent variable /. Its solution 
is U — a + 46/ + /_ (a, 6— real constants). Then functions 

Fi(x) = ki(a + at- exp(Ax) + /?+/?_ exp(— Ax)) 

+ k 2 (a 2 + a 2 _ exp(2Ax) + f3 2 + (3 2 _ exp(-2Ax)), 
— F 2 (x) = ki(a + f3~ exp(Ax) + exp(— Ax)) 

|| Due to Eq. (|215|) . the function F should be additionally representable in the form F — Fi(2xi) + 
F 2 (2x 2 ). 
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+ /c 2 (a+/r exp(2Ax) + (3j r a 1 _ exp(-2Ax)), 
c± = ± a ± exp(Ao;±) + (3± exp(-\x±) 
X(a± exp(Ax-i-) — f3± exp(— \x±)) 

(with k\ = aA 2 , k 2 = 46A 2 ) are real solutions of (I214p . ( I215p . if a±, (3± are real for 
A 2 > 0, and a± = (3* ± for A 2 < 0. 

5) To find a next class of solutions it is useful to rewrite (I214p in terms of Xi )2l 

(F 1 (2x 1 ) + F 2 {2x 2 ))d 1 {C+ + C_) + F{(2 Xl )(C + + C_) + F 2 '(2x 2 )(C+ - C_) = 0. 

Among known particular solutions the most compact one is, 

6 2 

= C-{x) = ax 2 + c, F x {2x x ) = 0, F 2 {2x 2 ) = — . (227) 

x 2 

After inserting these solutions (1222 j) - (12271) into the general formulas (12171) . one 
obtains the analytical expressions for potentials. Their explicit form can be found in 



Two additional classes of particular solutions of the system ( I203j) - (12051) were 
obtained analogously by means of suitable ansatzes for the case of degenerate metrics 
gik = diag(l, 0) ( see [79]) and the case of deformed hyperbolic metrics = diag(l, —a 2 ) 
with a^0 (see |86j). 

8.3. Integrability 

Similarly to the case of one-dimensional Polynomial SUSY Quantum Mechanics (Section 
3), partner two-dimensional Hamiltonians and can be considered as components 
of the diagonal super-Hamiltonian H, which are intertwined by components q* 1 of the 
off-diagonal supercharges Q ± . Thus, the operators H, Q ± constitute the deformed SUSY 
algebra: both supercharges are nilpotent operators, and they commute with the super- 
Hamiltonian. The difference is in the third relation of SUSY algebra: while in one- 
dimensional case anticommutator of Q ± produced the second order polynomial of H, in 
two-dimensional situation it gives some fourth order diagonal operator [76], [77], [78] . 

«™ -4 *-(T^)-(*r rV)" (228) 

In general, R(°\ are not expressed in terms of h^°\ in this case, but nevertheless 
they are related. Indeed, as it follows directly from the intertwining relations (I196p . they 
are the fourth order symmetry operators of corresponding systems, 

[hP\ RW) = [h {1 \ R<V] = 0. (229) 

In general, these symmetry operators might be expressed as functions of hS l \ i = 0,1, 
similarly to the one-dimensional case. But investigation [77], [79], [78] shows that this 
is not the case: two opportunities can be realized depending on the metrics g^ in the 
supercharges. 
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1) For the case of elliptic (Laplacian) metrics = 5ik (and only for this one), the 
fourth order operators R^\ i = 0, 1 are expressed in terms of Hamiltonians and the 
second order operators R?\ which also commute with hS l \ 

R (i) = (/ i »)2 + ^(i) + ^(i)^ = g (230) 

with n = Const. The explicit forms of differential operators R® and of corresponding 
Hamiltonians depend on the values of constants in (I206p . and they are given in 
|77j . [75] . It was shown that these Hamiltonians allow the standard procedure of 
R— separation of variables [143] in polar, parabolic or elliptic coordinates. 

2) For all other allowed metrics (I198p . the fourth order operators R^\ i = 0,1 
can not be reduced to any symmetry operators of lower order. This was shown in the 
general form in [77], and it follows also from the results of L.P.Eisenhart |146j . where 
the exhaustive list of systems which are amenable to separation of variables was given 
(just in Cartesian, polar, elliptic and parabolic coordinates). 

Thus, for all possible metrics ( 11981) each member of the variety of two-dimensional 
Hamiltonians h®, which satisfy the intertwining relations fl!96j) . has the symmetry 
operator and therefore is completely integrable. In general, this fact does not 
provide the solvability of the models, but some of them are partially or completely 
solvable (see the next Section). 

In papers [78], [79], the classical limit of Polynomial two-dimensional SUSY QM 
was considered. In particular, the prescription how to look for the integrals of motion 
(classical analogues of symmetry operators) of fourth order in momenta for given 
classical Hamiltonian was formulated. Namely, it is necessary to find such the complex 
classical functions q^(x,p) = {q^{x,p))^ in the phase space, that they satisfy the 
following relations: 

to h d } = if(x)q+; {q^, h d } = -if(x)q d , (231) 

where {■, ■} are Poisson brackets, and f{x) — an arbitrary real function of x. If such real 
function exists, the classical Hamiltonian h c i obeys the factorizable integral of motion 
Id = q^Qd of fourth order in momenta, 

{h d ,I d } = 0. (232) 
9. SUSY separation of variables for d — 2 

The idea of shape invariance was very productive in one-dimensional SUSY Quantum 
Mechanics (see Section 6 and references therein). In the present Section, the two- 
dimensional analogue of that shape invariance will be exploited effectively [80], [82J, 
|13j . [83] . [88] . [16] . [90] . [9~T] . [92] . [93] for two-dimensional Hamiltonians which satisfy 
the supersymmetrical intertwining relations (I196p . In contrast to one-dimensional case, 
there are essential differences in d = 2: in general, the ground state of h^(a) does not 
coincide with zero modes of second order supercharge q + (a), the value E (a) ^ 0, and 
also many zero modes of q + exist. Up to these crucial differences, all other actions can 
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be repeated. It was shown [80], [82], [13], [93] that such direct generalization of the 
procedure above may provide only partial (quasi-exact) solvability of the model (see the 
Subsection 5.2). But for some models [88], [90], [H], [92J it will help to solve the model 
completely (see Subsection 5.3). 

The idea to consider the quasi-exactly-solvable models - the intermediate class 
between exactly solvable and unsolvable analytically models - was introduced on 80- 
ties in [147] |153j . In particular, the series of papers by A.Turbiner, A.Ushveridze 
and M.Shifman was devoted to the elegant algebraic method of construction of 
one-dimensional quasi-exactly-solvable (and sometimes, of exactly solvable) quantum 
models. In general, this method works beyond the supersymmetry, but both approaches 
can be combined in one-dimensional case as in |154j . This approach is applicable to two- 
dimensional problems as well, but only in curved spaces with nontrivial metrics [151] . 

The realization of such opportunities to solve non-trivial two-dimensional models 
partially, or even completely, would be very important since only one regular method 
to solve analytically the Schrodinger equation for two-dimensional models is known: 
reduction to a pair of one-dimensional problems by means of the procedure of separation 
of variables |143] . This method can be used for very restrictive class of models, and full 
classification of models which allowed separation of variables was given by L.P.Eisenhart 
[146] : four possibilities exist - Cartesian, polar, elliptic and parabolic coordinates. The 
general form of potentials amenable to separation of variables is also known explicitly up 
to arbitrary functions of one variable. And analytical solution is possible only if these 
functions belong to the list of exactly solvable one-dimensional potentials. All these 
Hamiltonians H are integrable: the symmetry operator R of second order in derivatives 
(in momenta) exists, [H, R] = 0. Besides models amenable to separation of variables, 
the class of so called Calogero-like models [155] , [T8j is known as well. They describe 
the specific forms of pairwise interaction of N particles on a line, and they are solvable 
by means of special transformation of variables which leads to a separation of variables. 
Recently, new classes of solvable two-dimensional models were built in [156] . [157] . [158] . 
but all these models are superintegrable and amenable to separation of variables. 

In two following Subsections, we shall present two special procedures of SUSY- 
separation of variables in two-dimensional models, which do not allow to use the 
conventional separation of variables. Subsection 9.2 presents the first procedure of 
SUSY separation of variables where variables are separated in the supercharge. It leads 
to QES models, and the specific model of two-dimensional Morse potential illustrates 
this method. In Subsection 9.3 the second procedure of SUSY separation of variables is 
given where variables are separated in one of partner Hamiltonians. In the case of the 
same Morse model, but with particular values of parameter, it allows to solve the model 
completely, i.e. to find analytically the whole spectrum and all wave functions. 
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9.1. SUSY-separation of variables I: QES models 

The first variant of SUSY-separation of variables is realized when the Hamiltonians 
hP\ h^ 1 ' do not allow standard separation of variables, but the supercharge q + does 
allow (SO], H3], |85], |5Z!, US], [S3]. The general scheme is the following. Let's 
suppose that we know zero modes of q + , 

q + Q n (x) = 0; n = 0,l,...,N; q + Q(x) = 0. 

Due to intertwining relations f ll96p . the Hamiltonian h^ 1 ' obey the important property: 
the space of zero modes f2 n is closed under the action of 

h (1) Q{x) = CQ{x) 

with constant matrix C. If this matrix is known, and if it can be diagonalized, 

BC = AB; A = diag(X , Xi, X N ), 

the eigenvalues of h^ 1 ' can be found algebraically, 

h (1) {M(x)) = A(M(x)). 

Thus, to realize this scheme of construction some number of energy values and 
corresponding wave functions, we need: 

- to find zero modes Q n (x); 

- to find constant matrix B, such that BC = AB. 

As for the first step, the zero modes can be obtained by using the special similarity 
transformation (not unitary!), which removes the terms linear in derivatives from q + , 

q+ = e -^Ve +x{a ° = d\ - d\ + \{F x (2x x ) + F 2 {2x 2 )); 

X{x) = ~\{ J C + (x + )dx + + J C_(xJ)dx^j. 

Now, q + allows separation of variables for arbitrary solution of intertwining relations, 
and we obtain the first variant of new procedure - SUSY-separation of variables. 
Similarly to the conventional separation of variables, separation of variables in the 
operator q + itself does not guarantee solvability of the problem. 
The next task is to solve two one-dimensional problems, 

(-<9i - -F 1 {2x 1 )) , q n {xi) = e n 7? n (£i); 

{-dl + ^F 2 {2x 2 ))pn{x2) = e n p n (x 2 ). 
Three remarks are appropriate now. 

Remark 1: the same similarity transformation of does not lead to operator 
amenable to separation of variables. 
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Remark 2: the normalizability of Q n has to be studied attentively due to non- 
unitarity of the similarity transformation. 

Remark 3: we have no reasons to expect exact solvability of the model, but quasi- 
exact-solvability can be predicted. 

As for the matrix B, it must be found by some specific procedure. Such procedure 
was used in example which will be presented below. 

In principle, the first scheme of SUSY— separation of variables can be used for 
arbitrary models satisfying intertwining relations by supercharges with Lorentz metrics. 
The list of solutions of intertwining relations is already rather long (see Subsection 
8.2.2), and it may increase in future. The main obstacle is analytical solvability of 
one-dimensional equations, obtained after separation of variables in the operator q + . 

Below we describe briefly such a model which can be considered as the generalized 
two-dimensional Morse potential, 

(XX— 

C+ = Acta: C- = Aaa ■ coth ■ 



2 ' 

fiixt) = \F % {2 Xi ) = -A(e- 2ax >-2e- ax >); i = 1,2 



0°)>« = a 2 a(2a T 1) sinh" 2 (^) + 4a 



V 



+ A 



e ~2axx _ 2e~ aXl -|- e ~ 2ax 2 _ 2e~ aX2 



(233) 



where A > 0, a > 0, a— real. 

To explain the name, we present the potentials in the form, 

V(x) = V Morse ( Xl ) + V Morse (x 2 ) + v(xi,x 2 ), 

where first two terms are just one-dimensional Morse potentials, and the last term mixes 
variables xi, x%- 

The solutions of one-dimensional Schrodinger equations with Morse potentials 
VMorse(x) are well known [159] , and the zero modes can be written [80"] . [T3] as, 

a n (x) = (^.^2_) 2a exp (-ii±^)(e 1 e 2 ) s -F(-n, 2s n +l- e x )F(-n, 2a B +l; (234) 
v V A |?2 — ?i| 7 1 

c 2 v^4 , , VA 1 

£i = exp(-axi); s n = n - - > 0. (235) 

a a 2 

The conditions of normalizability and of absence of the "fall to the center" are, 



11, VA 1 

a 6 (-oo, - ^=); Sn = — -n-->-2a>0 

To obtain the matrix C explicitly, one must act by on Q n . The matrix turns 
out to be triangular, and therefore, the energy eigenvalues coincide with its diagonal 
elements, 

E k = c k k = -2(2aa 2 s k - e k ). 
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To find a variety of wave functions is a more difficult task. For that it is necessary 
to find all elements of C and all elements of matrix B. The recurrent procedure for the 
case of two-dimensional Morse potential was given in [80], [T5] . This variety can be 
enlarged by means of shape invariance property (1127ft of the model, which can be easily 
checked. Similarly to one-dimensional shape invariance, each wave function constructed 
by SUSY-separation of variables leads to a set of additional wave functions, which can 
be written as, 

q (a)q (a--)...q (a — )*(o - — 



frw M ^ M-l. 

[E (a - —) + TZ{a — ) + ... + U{a 



x 



q (a)q (a--)...q {a — )*(a - — 7J 



with integer M. Analogous approach works for the two-dimensional generalization of 
Poschl- Teller model [82], for some two-dimensional periodic potentials [87] and for the 
two-dimensional generalization of Scarf II potential [93] . 

9.2. SUSY-separation of variables II: Exact solvability. 

Among all known solutions of two-dimensional intertwining relations with second order 
supercharges a subclass exists [88], [90], [M], [16], [92], where one of intertwined 
Hamiltonians is amenable to standard separation of variables due to specific choice 
of parameters of the model. Its superpartner still does not allow conventional 
separation of variables. 

The scheme will be described below for the same specific model - two-dimensional 
generalization of Morse potential, 



V®>W = a 2 a(2a T 1) sinh" 2 (^) + 4a V + A 



g— 2axi 2e~ axi + e~ 2aZ2 2e~ 01x2 



2 

Let's choose ao = —1/2 in order to vanish the mixed term in and therefore, 

the Hamiltonian allows the conventional separation of variables. Moreover, after 
separation of variables each of obtained one-dimensional problems is again exactly 
solvable: they coincide with one-dimensional problems of the previous Subsection, where 
they occured in a different context (separation of variables in q). 
The discrete spectrum of this one-dimensional model is, 

2 2 VA 1 

e n = -a s nl s n = n-->0; n = 0,1,2,..., 

a Z 

Wave functions are expressed in terms of degenerate hypergeometric functions, 



£j 2v A 

T] n (xi) = exp(-^)(&) s "F(-n, 2s n + 1; & = exp(-aXi). 

2 a 

Due to separation of variables, the two-dimensional problem with h^\x) is exactly 
solvable. Its energy eigenvalues are, 

E n m E m n £ n -|- 6 m , 
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being two-fold degenerate for n 7^ m. The corresponding eigenfunctions can be chosen 
as symmetric or (for n/m) antisymmetric combinations, 

^fiiff 4 ^) = Vn(xi)r] m (x 2 ) ±rj m (x 1 )r] n (x 2 ). 

Our initial aim here is to solve completely the problem for h^°'(x) with ao = —1/2. 
The main tool is again the SUSY intertwining relations, i.e. isospectrality of and 

but up to zero modes and singular properties of q ± . In general, we may expect three 
kinds of levels of h^°\x), 

(i) . The levels, which coincide with E nm . Their wave functions can be obtained 
from by means of q + . 

(ii) . The levels, which were absent in the spectrum of h^(x), if some wave functions 
of h(°\x) are simultaneously the zero modes of q~ . Then the second intertwining relation 
would not give any partner state among bound states of h^(x). 

(iii) . The levels, which were also absent in the spectrum of h^(x), if some wave 
functions of h^°'(x) become nonnormalizable after action of operator q~ . 

We have to analyze these three classes of possible bound states of one after 
another. 

(i) . The first SUSY intertwining relation gives the two-fold degenerate wave 
functions of with energies E nm , ty^ = g + \l/^ m . But q + includes singularity on the 
line x\ = X2, therefore the normalizability of depends crucially on the behavior 

of ^E nm 011 the line £1 = £ 2 . One can check that only antisymmetric functions yfr^ 
survive, i.e. only symmetric survive. This fact can be demonstrated [88] both by 
direct calculation and by indirect method - by means of symmetry operator R^°\ 

The indirect method explores that the symmetry operator = q~q + for 

ao = —1/2 can be written in terms of one- dimensional Morse parts hi(x\), h 2 (x 2 ) of 
the Hamiltonian = h\(xi) + h 2 (x 2 ), 

r(o) = (h^xx) - h 2 (x 2 yj + 2a 2 (h 1 (x 1 ) + h 2 (x 2 yj + a 4 . 

Therefore, 

R {0) V { ZiW = rn, m y£ n Jx); r„, m = a%n - mf - l][(s n + s m f - 1], 

and 

II 11 2 = im^ A 1 a -q + 1 ¥° )A ) = r \\¥° )A 11 2 

II Ej n ,m II > E n , m I ^ ^ I c n , m .' n,m\\ E n . m \l 

For n = m, wave functions ^E~!f„ van i sn identically by trivial reasons. It is clear 
now that wave functions ty^ S also vanish. For all other n,m, functions ty^ S have 
positive and finite norm, and there is no degeneracy of these levels. 

(ii) . These possible bound states of are the normalizable zero modes of q~. 
The variety of such zero modes is known from [80] : they exist only for positive values 
of a x 1 

^^ + ^,+00), 
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which does not contain the value ao = —1/2. Thus, no normalizable bound states of this 
class exist for h(°\ 

(iii). We have to study an opportunity that q~ destroys normalizability of some 
eigenfunctions of h^°\ It could occur due to singular character of q~ at x\ = x%. The 
analysis was performed [88] in suitable coordinates. It shows that q~ is not able to 
transform normalizable wave function to nonnormalizable. Therefore, the third class of 
possible wave functions does not exist too. 

Summing up, the spectrum of with a = —1/2 consists only of the bound 
states with energies E nm for \n — m\ > 1. This spectrum is bounded from above by 
the condition of positivity of s n , s m : n,m < yA/a — 1/2. The corresponding wave 
functions are obtained analytically |88j . 

The results above can be expanded to the whole hierarchy of Morse potentials 
with a k = —(k + l)/2 with k = 0,1, ... by means of shape invariance property. Let's 
denote elements of the hierarchy as h^(x;a k ), h^(x;a k ). All these Hamiltonians are 
also exactly solvable due to shape invariance of the model, 

h { -°\x-a k ^ l ) = h {1 \x-a k )- k — l,2, 

This means that the following chain (hierarchy) of Hamiltonians can be built, 

h m {x; a )^°\x; a ) = h {1) {x; ai )+h m {x; oi) = . . .+h {1) {x; a k - X ) = h^{x; a k )^h w {x; a k ), (236) 

where the sign denotes intertwining by q ± (a,i). 
In the general case, the functions 

^S, m ( f ' a fc) = ? + ( a fc)*£, m (^ a fc) = g + (a fe )g + (a fe _i) . . . g + (a )^^(x; a ) (237) 
(if normalizable) are the wave functions of h^(x; a k ) with energies E nyTn = — a 2 (s^ + s^). 
The symmetries of wave functions alternate and depend on the length of chain. This is 
true but up to zero modes of operators q + . 

It is necessary to keep under the control normalizability of \l/ and zero modes of q + . 
This control is performed algebraically by means of identity, which must be fulfilled up 
to a function of H, 

RW(a k ) = R<-°\a k - 1 ). 
Actually, the following equation holds, 

q-{a k )q + {a k ) = q + (a k ^)q~ (a fc _i) + a 2 (2k + 1) ^2h {0 \x; a fc _i) + a 2 (2k 2 +2k + 1)] .(238) 

These relations allowed to evaluate the norms of wave functions. The result is the 
following. The spectra of Hamiltonians h^(x;a k ) are not degenerate. They consist 
of the bound states with energies E n>m , with indices \n — m\ > k + 2, and their wave 
functions (x;a k ) were given analytically above. 

The procedure of this Subsection was reproduced fully for two other two- 
dimensional models which are not amenable to standard separation: generalized Poschl- 
Teller [90] and Scarf II [92] potentials. The full energy spectra and corresponding 
wave functions were also built analytically by the second variant of SUSY-separation of 
variables. 
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10. Retrospectives, perspectives and missing points 

The purpose of this review has been to elucidate the recent progress in Nonlinear 
SUSY realization by intertwining Darboux transformations for a community working 
on methods of isospectral design and their applications. We wanted also to draw an 
attention to new (quasi-)exactly solvable potential systems in two dimensions non- 
separable by a conventional choice of coordinates. Summarizing our experience we 
describe the general SUSY QM as governed by the extended nonlinear SUSY algebra 
with J\T pairs of nilpotent supercharges Qj, Qj = Qj = and a number of hermitian 
hidden-symmetry differential operators R a = R^, [R a , Rp] = 0; < a, f3 < M. 

Not all the research areas that are linked to the subject of the present review 
have been covered in the main text. On the other hand, there remains a plenty of 
open questions and challenges to be solved. Below we list some of these directions and 
mention few missing topics. 

• The method of multidimensional supersymmetric intertwining relations and shape 
invariance turned out to be useful for study of different variants of Calogero-like 
models of N particles on a line, including the models with internal degrees of 
freedom and with pairwise interactions based on any root systems (95] , [160] - |167] 
(see the recent review |18j). 

• One-dimensional supersymmetric models with matrix superpotentials and 
supercharges |168j |173] have been used either for description of motion of 
spin particles in external fields [65] [67], [97], jlOlj . |174] [179] or to study 
scattering of particles with strong coupling of channels [174] . [175] . |180j [183] 
with practical application of the results to specific processes. Such systems are not 
fully investigated in searching of extended SUSY systems with hidden symmetries. 
It is clear that a comprehensive understanding of irreducible building blocks for 
spectral design analogously to the scalar case is not yet achieved. Also, matrix 
models with higher dimensionality of space have to be studied. To start with, 
two-dimensional models with matrix superpotential were studied in [184] . 

• Coherent states in the framework of SUSY QM models were built in [185] - [188], 
but only for one- dimensional systems [189]. 

• Formulation of multidimensional SUSY QM in arbitrary curvilinear coordinates 
was developed in [190] (see also [191] - [193] ). This formalism could be useful for 
investigation of SUSY design of different cosmological and brain world models. 

• Generalization of Polynomial SUSY QM for higher (first of all, d = 3) space 
dimensions seems to be very interesting (see [89]). 

• Effective mass SUSY QM is a rather popular generalization of conventional SUSY 
QM for the Schrodinger equation with coordinate-dependent mass. This sort of 
models may have various applications in physics of solid state and nano devices as 
well as in cosmology. The well developed SUSY approach to effective mass QM 
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is rather promising for solution of different problems in physics. Some of suitable 
references are |194] - [200J. 

• The periodic potentials were studied actively in the framework of SUSY QM as 
well, see |2Ulj - |2U8] , jSH] on one-dimensional models and |87J for two-dimensional 
quasi-exactly-solvable models with periodic potentials. Certainly, spectral design of 
periodic potential systems is promising, especially, in two-dimensional QM related 
to new nanoscale materials like graphene. 

• The embedding of shape invariant potentials and their eigenfunctions into so called 
Discrete Quantum Mechanics has been analyzed recently [209] . (210J. 

• The generalizations of conventional SUSY QM to the models with complex 
potentials [HI] , [HI] , [211] - [213] provide an effective tool in study of non-Hermitian 
[214J [219] and, especially, PT symmetric [220] [225J Quantum Mechanics 
with real energy spectrum (see reviews [226] . [227] ). Spectral design of non- 
Hermitian systems with non-diagonalizable Hamiltonians was investigated both in 
one- dimensional |228[ \229\ 1230] and two-dimensional [231] cases. 

• The polynomial SUSY in two dimensions has already brought a number of examples 
of new type of irreducible SUSY with hidden symmetries of higher-order in 
derivatives [76], [77]. One may expect a variety of new types of irreducible SUSY 
for third-order (and higher-order) supercharges as well as new discoveries in two 
and three dimensions. 

• The interesting direction for investigations is the interrelations between the inverse 
scattering problem and methods of SUSY QM [232] - [237] . 

• The SUSY QM approach happens to be applicable and useful for building of 
different brane-world models with extra space dimension [238J - [241J. 
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